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1.0 INTRODUCTION 


The objective of this report is to describe a computer generated 
false color display program which creates digital multicolor graphics. 
The primary goal of this program is to display geophysical surface 
parameters measured by the SEASAT-A satellite scatterometer (SASS). 

This False Color technique shall be used to assist the NASA personnel at 
Langley Research Center involved in the SASS experiment. 

The SEASAT-A satellite has a special purpose microwave radar scat- 
terometer which is used to measure the normal ized radar cross section 
(NRCS) of the earth's surface. While operating in its various modes, 
the satellite can measure NRCS for fifteen angles of incidence both 
forward and aft of the flight path. The cross-section measurements are 
used to infer surface roughness and reflectivity. Over the ocean the 
NRCS can be used to calculate surface wind speeds and define ice/water 
boundaries [ref.l]. Overland it has been conjectured that the NRCS 
could demonstrate seasonal variations in soil moisture [ref. 2]. 

False color displaying of data sets is simular in principle to con- 
touring or gray-scal ing. The data is incrementally scaled over the 
range of acceptable values and each increment and its data points are 
assigned a color. The advantage of the false color display is that it 
visually infers cool or weak data versus hot or intense data by using 
the rainbow of colors. For example, with wind speeds, levels of yellow 
and red could be used to imply high winds while green and blue would 
imply calmer air. 

SASS data is sorted into geographic regions and the final false 
color images are projected onto various world maps with superimposed 
land/water boundaries. FiQures 1.1 and 1.2 are two examples of the SASS 
data taken during orbit 331, and are shown here in a polar stereographic 
and mecator map format, respectively. It is this final format that will 
provide an excellent qualitative as well as quantitative look at the 
SASS sensor measurements. 
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Figure 1.1 The final display format using a Polar Stereographic map projection. 
Data is from beams 2 and 3 of REV 331. 
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Figure 1.2 The final display format using a Mercator map projection. Data 

is taken about Hawaii during Hurricane Fico from beam 2 of REV 331. 
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1.1 Requirements 

SEASAT-A data is received and processed by Jet Propulsion Labora- 
tories which then issues the appropriate sensor data on mag tapes to 
various contractors. At the time of this report, a complete tape with 
the final format was not available. However, several tapes have been 
issued which have partial but adequate sensor data records for 
demonstration of the program 

It was a requirement that the final program be adaptive in the 
sense that the user should be able to control data selection and the 
final false color format. The basic specifications imposed on the 
software can be listed as follows: 

• The SASS sensor records shall be read and sorted by mode, 
beam number, polarization and cell number into earth- 
located data subsets. 

• The data subset shall be gridded using selectable pixel 
sizes, and an optional pixel blending feature shall be 
used to interpolate for missing SASS data. 

• The gridded data shall be displayed in either a mecator or 
polar stereographic map format. 

• The display color scale shall have a selectable linear, 
logarithmic or exponential relationship to the data. 

• The final display shall be superimposed over land/water 
boundaries. 

1.2 Outline 

The handling of the SASS data and the ultimate creation of a false 
color display is handled in a self-contained routine called FALSCLR. 
Called by the main routine are several subroutines which use display 
algorithms to create the false color data file. The purpose of Section 
2.0 is to describe the various gridding, smoothing and map projection 
al gorithms. 

It became apparent during the course of the software development 
that four separate programs were necessary to fulfill the task require- 
ments. Besides FALSCLR, two routines were developed to handle display 



annotation and land/water boundaries, respectively. A fourth program 
was developed to merge these three display files into a complete false 
color picture. Section 3.0 is devoted to the description of each of 
these four routines. 

Section 4.0 contains a description of how to coordinate the four 
routines into a meaningful job. The section contains a complete outline 
of program operation so that the routines can be run by the user who is 
reasonably familiar with the false color equipment. Since job control 
cards are an important part of program control; and since it is possible 
to run multiple programs with the proper cards, a discussion of 
candidate control cards is included. Also, described in detail are all 
possible data cards along with their input formats. 

A discussion on the operation of the digital display system is 
contained in section 5.0. Included is a brief description of the 
display equipment and a simple guide to its operation as it might apply 
to the false color display techniques mentioned here. The user may find 
this section useful in his initiation to the system. 

Section 6.0 is devoted to a discussion of future development in 
program software. It will be necessary to eventually adapt the routines 
to handle new data tapes as they arrive. It has also been suggested 
that a routine should be developed to do selective normalization of the 
SASS data. These and future modes of operation are the subject of this 
section. 
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2.0 FALSE COLOR ALGORITHMS 

The purpose of this section is to detail the algorithms used by the 
false color program to generate packed integer map projection files from 
the raw SASS data. All the following algorithms are executed within the 
framework of the basic raw data to display routine called FALSCLR. The 
algorithms are very general and can easily be adapted to similar 
problems where data must be digitized. 

2.1 Data Gridding and Blending 

The SASS measurement of the normalized radar cross section (NRCS) 
is taken along swaths of area over the earth that correspond to about 
950 kilometers along either side of the satellite ground flight path. 

The data is taken continuously as the satellite orbits around the earth, 
and the final data tapes contain segments from these orbits. Once the 
data has been sorted into geographic subsets a large file of unequally 
spaced data points remains. To digitize these SASS data, one first 
extrapolates equally spaced points within the desired geographic region. 

Then the proper SASS data value is assigned to each of these points by 
interpolation. The operation is referred to as data gridding and is 

performed in a routine called SSGRID, 

The gridding routine is designed to project unequally spaced points 
onto a grid matrix which is a linear 1 atitude-and-longitude matrix of 
equally spaced grid boxes called pixels. The pixel size in kilometers 
is a variable controlled by the user. The final grid matrix has dimen- 
sions n X m given by: 

n = * 111.112 km/pixel size 

i ^ (2-11 

, 1 , = * 111.112 km/pixel size ' ' 

where ^2 <^efine the latitude and longitude coordi- 

nates, respectively, for the region of interest, The angular size of 
each pixel is defined as: 

A^= (’^1- 

(*i - *2)/n (2-2) 
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Each pixel is assigned a data value by computing the weighted average 
of data points within the pixel. The weighting used is known as inverse 
squared distance or gravity weighting. Only those data values which 
fall within a particular pixel can contribute to the value of that 
pixel. 

The SASS data file contains unequally spaced data points U|^(a|^, 

'^ 1 ^. where a|^, are the SASS parameter value, lati- 

tude position, and longitude position, respectively, for each data 
point k. The gridding algorithm is defined as follows: 

For each data value U|^ in the SASS data file, perform the 
following functions: 

a) compute the matrix address (i,j) of the grid pixel containing 
this point by: 

i = 1 

( 2 - 3 ) 

j = 1 


b) define the coordinates of the grid pixel (i,j) as x^-j, 
y^-j and compute the dimensionless range r|, as: 




( 2 - 4 ) 


c) compute the running sums P(i,j) and Q(i,j) as follows: 


Q(i,j) =5Z l/'^k 
k 

After the entire data set has been exhausted, the pixel data values 
V(i,j) are given by: 


V(i J) = P(i,j)/Q(iJ) 


(2-6) 
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If no data falls within a particular pixel, that pixel is flagged for miss- 
ing data. 

Associated with the gridding process is an optional pixel blending 
feature. Upon request, a routine called SMTHGD is called to interpolate 
data values for pixels which were flagged for missing data. The pixel value 
is computed as the weighted average of the pixels surrounding that particu- 
lar pixel. As in the gridding routine, gravity weighting is used to inter- 
polate this value. The degree of blending is user controlled. The value of 
the grid pixel (i,j) is found to be: 

V(i,j) = 2 . V(k,ji) (2-7) 

tk 


where (k,£) is the 2-dimensional address of the pixels, (1) surrounding the 
pixel (i,j) and (2) having range r|^^ from the center of this pixel. 

are the normalized weights associated with the pixels V(k,£) 
having range r|^^. 


2.2 Map Projections 

The map projection algorithms are designed to perform the non-linear 
transformation of gridded SASS data onto two common map types. The mercator 
and polar stereographic maps were selected because they complement each 
other in areas of minimum land mass distortion, and SASS data analysis can 
be effectively done with these map types. 

2.2.1 Mercator .- The mercator map is a cylindrical projection of the 
earth's surface and is ideally suited for the tropic latitudes [ref. 3]. 

Full world projections are possible with this map type; however, consider- 
able vertical map distortion will result. The properties of the mercator 
map are that the longitude 1 ines are parallel and are of equal distance 
throughout the map. The latitude lines are also parallel, but their spacing 
is a function of latitude. The coordinates of the map are those used for 
the grid matrix, and the map projection is scaled to fit within the user's 
specified pixel dimensions. 

The algorithm for projecting the grid matrix onto the mercator map can 
be throught of as a magnifying routine which stretches or shrinks the grid 
matrix to fit the map. The longitudinal or grid column magnifying factor Ay 
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is constant throughout the transformation and is given by: 


Ay = [nyn] 


(2-8) 


where n^ and n are the longitudinal pixel sizes of the map and grid, 
respectively. The latitudal or grid row magnifying factor is a loga- 
rithmic function of latitude and is calculated for each row in the grid 
matrix. First, a scaling factor and the map's projected equator 
distance are given by: 


s = mj/(a-b) 
eq = m 2 *a/(a-b) 


(2-9) 


where m^ is the map's latitudal pixel size and a and b are given as 

a = logg(tan(45° + *^^/2)) 
b -- logg(tan(45° + '/'2/2)) 


( 2 - 10 ) 


Then, for the jth row of the grid matrix, we find: 
Cj = logg(tan(45° + - j.Ay)/2)) 


( 2 - 11 ) 


The jth latitudal magnifying factor is given as: 


AXj = [eq - dj] 


( 2 - 12 ) 


The map values M(i,j) are finally given as: 


For k=l, ... , AXj and j = 1, ..., Ay 


(2-13) 


M(k,il) = V(i,j) 


V'/here V(i,j) are the ith, jth grid values. 
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2.2.2 Polar stereographic .- The polar stereographic map is a plane 
conformal-azimuthal projection that is ideally suited for the polar 
regions. Because it is a continuous map of the globe, it is more useful 
in hemisphere charts than the mercator. The polar stereographic map has 
latitude lines which are closed concentric circles about the pole and 
the longitude lines which radiate outward from the center. The map is 
scaled to fit within user inputted dimensions and has an outermost 
latitude ’^3 such that *^3 ± 

Projecting the grid onto the map is completely non-linear and is 
considerably more complicated than the mercator. Basically, the algor- 
ithm calls for calculating the center latitude and longitude of every 
pixel in the map matrix thereby determining which grid pixel is associ- 
ated with it. The calculation goes as follows: given an x \2 x r \2 map 
matrix calculate the following scale factor: 

sj = l/ 2 -n 2 /tan( 45 ° - ’^ 3 / 2 ) (2-14) 

where ^3 is the standard or outermost latitude. For each Ath, kth map 
pixel compute 

r2 =(03 - + (h 3 - k)2 (2-15) 


where 03 = 02/2 and r is the radial distance to the pixel from the 
map's center. Next compute the pixel's latitude and longitude 
as follows: 

= 2*tan"^(r/Sj) (2-16) 


= tan“l((n3 - k)/(n3 - ii)) + 180° 

The grid matrix element address (i,j) which corresponds to the ,Uh, kth 
map element is found as: 

i = (^ - 

(2-17) 

j = C'-i - 
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Finally the map value M(ii,k) is given as 


= V(i,j) (2-18) 

where V(i,j) is the ith, jth grid value. 

2.3 Display File Preparation 

To display the map projection on the digital display equipment, a 
special raster file must be created on magnetic tape. This file will 
contain packed integer values which have a functional relationship to 
the map's SASS parameter values. A color scaling algorithm is designed 
to perform the transformation of the map values into scaled integer 
values which ultimately represent different colors. 

It was decided that at most 20 unique colors should be used at any 
time in the false color display. These colors will have either a 
linear, logarithmic or exponential relationship to the SASS data. The 
color algorithm goes as follows: 

Define an upper bound M-| and lower bound M2 map data 

values M(i,j), and let be the number of colors desired. Compute 
AM, the color incrementing value in map data units as follows: 

Linear scale: 


am = (Mj - M2)/Ne 

Logarithmic scale: 


(2-19) 


am = N(,/Togg (Mj - M2 + 1) 
Exponential scale: 


( 2 - 20 ) 


am - (M;^ - M2)/logg (N^ + 1) 


(2-21) 
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For each ith, jth map matrix element compute a color integer value 

I 

N(i ,j) as follows 


If M(i ,j) > 

If M(i,j) < M 2 , 


1 et N(i J) = 

1 et N(i ,j) = 1 


( 2 - 22 ) 


If M(i,j) is within M^ and M 2 then, depending on the functional re- 
lationship, N(i,j) is chosen as follows: 


Linear scale 

N(i ,j) = (M(i ,j) - M2)/AM + 2 


(2-23) 


Logarithmic scale 

N(i,j) = AM*logg(M(i ,j) - M 2 + 1) + 2 


(2-24) 


Exponential scale 

N(i,j) 1 + e (M(i,j)-M2)/AM (2_25) 


For map values which were flagged for missing data, zero is assigned to 
N(i,j). 

Once the integer assignment is complete, it is necessary to pack 
the array N(i,j) such that only eight bits are used to represent each 
integer on the magnetic tape file. This eight bit form is necessary 
since the tape drive is controlled by a Kennedy model 9218 format 
control unit which is IBM compatible. Since the CDC computer uses a 
standard 60 bit word, it is necessary to pack 7 1/2 eight bit words into 
each word in the output list. This is accomplished by obtaining the 
lower eight bits of each N(i,j); bit shifting each word a multiple of 
eight bits; and packing these final, reconstructed words into an output 
array. This final array is written to tape using free format to 
guarantee that all 60 bits are transferred. 
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3.0 PROGRAM DESCRIPTION 


This section is devoted to the description of the entire false 
color display software. Four distinct program units were developed to 
ease program execution and reduce time and cost by eliminating repeti- 
tive runs. The first three units each generate a display file 
containing part of the total picture. The fourth program unit merges or 
superimposes these partial images into the complete digital display. 
Figure 3.1 is the simplified flow diagram describing the total job. A 
listing of each of these routines is found in the Appendix. 

3.1 SASS Display Generation 

Routine FALSCLR is a Fortran program designed to process and graph- 
ically display various SASS data sets. The program has two logical 
inodes of operation which are defined here as preprocessing and postproc- 
essing. These tv/o modes can be executed jointly or separately at the 
users descretion. The program's final output is a digital display file 
containing integer-represented SASS values projected onto one of two map 
images. Refer to the job flow diagram in Figure 3.2. 

3.1.1 Preprocessing .- In the preprocessing mode of operation, 
FALSCLR accesses the SASS data tapes and performs data sorting and 
selecting. These SASS tapes contain segments of the satellite's revolu- 
tions; and within each segment, there are unique sensor data records. 

Each record contains the mode of operation defined by the mode, 
polarization and beam number, the satellite's geographic location, and 
the fifteen measurements of NRCS taken from the particular antenna in 
use. A subroutine called DATAIN was designed to read these SASS tapes 
and pick out the appropriate data. Upon execution, this subroutine 
prompts the user to define the geographic coordinates of the desired 
data subset, to enter the desired mode of operation and to specify which 
of the fifteen NRCS measurements are needed. After the tape search is 
complete, the final data subset is written to a working file v/hich is 
accessable to the postprocessing routines. 
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FALSCLR 


A 



Figure 3.1 Simplified flow diagram of the complete job. 
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Figure 3.2 Flow diagram of program FALSCLR. 
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3.1.2 Postprocessing ." The postprocessing mode contains routines 
to access the preprocessed data subset and create the final display 
file. There are six subroutines invoked by FALSCLR to perform this 
operation. The algorithms used by these routines are discussed in 
detail in the previous section. The final output file is format 
compatible with the digital display equipment and is stored on disk for 
future access. 

Two of the six subroutines, SSGRID and SMTHGD, are designed to 
access the data subset and create an equally-spaced grid matrix of SASS 
values. Subroutine SSGRID requires the user to define the geographic 
bounds of the grid and the size in kilometers of the grid's square 
pixels. The routine then develops a linear latitude-and-longitude grid 
matrix within the requested bounds. SASS values are then assigned to 
each grid pixel by weighted interpolation of all the SASS data 
geographically located within the pixel. 

At the option of the user, a grid smoothing routine SMTHGD is in- 
voked. This routine attempts to define a SASS value to those grid 
pixels containing no SASS data points. The interpolation is done by ob- 
taining the weighted average of the pixels surrounding that pixel. The 
smoothed grid matrix is written to a working file which is passed on to 
1 ater stages of the program. 

The next two subroutines in the postprocessing operation, MERCAT 
and POLAR, were designed to operate exclusively of one another. They 
both access the above grid matrix and project the data values onto a 
nonlinear map matrix. If invoked, subroutine MERCAT obtains the grid 
file and projects this file onto a log-linear matrix defined by the 
equations of a mere ator map projection. The mercator map is a cylindri- 
cal projection which has linearly spaced longitude lines and logarith- 
micly spaced latitude lines. The size of the map is user specified and 
should be such that minimal area distortion is maintained. 

Paralleling MERCAT is subroutine POLAR which projects the gridded 
data onto a map matrix defined by the equations of the polar 
stereographic map. This map type is a pi ane conformal projection which 
has concentric rings of latitude, about either pole, spaced tangentially 
from the map's center. The longitude lines are equally spaced and 


- 16 - 


radiate outward from the center. The map matrix dimensions are square 
and, though user specified, there should be some constant value for all 
projections of this type. 

The final two subroutines in the postprocessing mode are routines 
COLORS and PACKEM. These two programs create the proper display file 
for projection on the digital display equipment. Subroutine COLORS 
accesses the map file described above and redefines the SASS values as 
integers. These integers have either a linear, logarithmic or 
exponential relationship to the SASS data. There are at most twenty 
integer values which ultimately represent different colors in the false 
color displ ay. 

Following COLORS is subroutine PACKEM which creates an output file 
which is a packed version of the file created by COLORS. To be compat- 
ible with the digital display tape drive, the above integer values 
should be written to tape using only the lowest eight bits. Routine 
PACKEM accomplishes this by packing several CDC 60 bit words with 7 1/2 
eight bit integers. It then outputs the 60 bit words to magnetic tape. 
This tape file can be displayed on the digital CRT for test purposes; 
however, the picture is not complete until this display is matched with 
those from the other program units. 

3.2 Display Annotation 

Routine ANNOT is a Fortran program designed to create annotation 
for the final false color display. ANNOT generates a plot vector file 
which contains non-standard plot calls. This plot file is eventually 
converted to a special raster plot file which is compatible with the 
display equipment. This file is saved on disk to be merged later with 
the false color file generated in FALSCLR. Since many displays would 

likely use the same annotation, this display file could be used repeat- 

1 

edly at a substantial savings of time and cost. Refer to the flow dia- 
gram in Figure 3.3 

The annotation generation by ANNOT is controlled by the user at 
the time of execution. A ten character map code is entered which 
generally defines the final displays using this annotation. Other 
various control values are input which define the color legend. This 
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Figure 3.3 Flow diagram of routine ANNOT. 
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legend contains color rectangles defining the total band of colors used 
and the color boundary data values in decibels. There is also a map 
label which states that this is a map of normalized radar cross section 

and that the scale is either linear, logarithmic or exponential in 
decibels. 

The annotation plot file is generated by invoking various plot sub- 
routines from the Langley Research Center's (LRC) graphics package. 

These subroutines create a plot vector file that generally would be 
accessed by the CalComp plotter for plotting the file image onto paper. 
However, routine ANNOT uses several non-standard CalComp plot commands 
which make the file incompatible with the CalComp plotter. Instead of 
plotting the file, a special raster generation routine called RASTGEN is 
implemented which converts the plot vector file to a raster type file. 

This raster file is saved on disk for future merging with the other 
display files. 

Program RASTGEN was developed at LRC by Robert Oakes and revised 
for general usage by Frances T. Meissner. It is used to convert any 
plot vector file containing the proper non-standard CalComp plot calls 
into a device-compatible image file for displaying on the digitizing 
equipment. RASTGEN is located on disk and is available to anyone who 
knows the proper user number. This disk file is a source image and must 
be compiled before execution. The raster file created by the program 
always contains 512 records each having 70 words. This translates to 

512 X 512 optical pixels when seen on the CRT or photographed by the 
Dicorned imaging recorder. 


3.3 Map Outline Generation 

Routine MAPGEN is a Fortran program which generates a plot file 
containing map images. These map images are either polar stereographic 
or mercator and contain latitude and longitude lines and superimposed 
land/water boundaries. The plot file is converted to a special raster 
file as was done for program ANNOT. The raster file is saved on disk 
for future merging with the other display files. Refer to Figure 3.4 
for the job flow diagram. 
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COMPUTE MAP 
CENTER LAT. 
AND ANGULAR 
LIMITS 



CALL SUPMAP 
TO PLOT 
MERCATOR 
MAP 



Figure 3.4 Flow diagram of routine t'lAPGEN. 
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MAPGEN creates a plot file by making use of a special graphics 
package available at LRC. This package is called SUPMAP and was 
developed by the National Center for Atmospheric Research (NCAR) in 
Boulder, Colorado. The routines in this package are designed to create 
a plot file which is compatible with the Cal Comp plotter. The SUPMAP 
plotting options available to user are : (1) a geographic region can be 
plotted using one, of nine possible map projections; (2) if desired, 
land/water boundaries are plotted using dashed or solid lines, (3) a 
border can be drawn around the plot frame; (4) U.S. state borders can be 
plotted; and (5) the latitude and longitude lines can be plotted at 
selectable intervals. MAPGEN makes use of all but option 4. The lati- 
tude and longitude line intervals are chosen to be ten degrees but can 
easily be changed. The plotted land/water boundaries have a specified 
resolution of one degree which is the one limiting factor v/hen genera- 
ting detailed map projections. 

SUPMAP has been revised and adapted to the LRC facilities by 
Patricia A. Winters. A new plot call translator was developed which 
allows the user to scale the plot image and offset the plot origin. 

This option was necessary so that the final plot image is identical in 
scale and location to the data display file created in FALSCLR. The 
plot translator also includes various non-standard CalComp plot calls 
similar to those mentioned in ANNOT. 

The plot vector file created by MAPGEN is temporarily stored on 
disk. Execution of RASTGEN completes the necessary conversion of this 
file into a display image raster file, much the same as in ANNOT. This 
raster file also contains 512 records each having 70 words. In this 
form, it is compatible with the display equipment; and, if projected 
onto the CRT, uses 512 x 512 optical pixels. Once generated this map 
file coqld be used repeatedly for displays requiring the same map type 
and area. The file is saved in the user's disk area along with the 
files created by FALSCLR and ANNOT. Together they complete the total 
fal se color displ ay. 
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3.4 Merging 

Routine MERGE is a Fortran program designed to do a bit-^by-bit 
merge of two binary files. The routine is intended to superimpose two 
display files and output this result to a new file. If run twice using 
the files created by FALSCLR, ANNOT, and MAPGEN, routine MERGE can 
create the complete false color display of the SASS data. Refer to the 
flow diagram in Figure 3.5. 

To execute MERGE, two raster display files must be available to the 
program, i At least one file must have 512 records, such as those created 
by ANNOT and MAPGEN. If the second file is smaller than this, an offset 
is computed which properly aligns the eight bit integer values of the 
two files. This offset is necessary for all files created by FALSCLR, 
which are at most 400 records by 54 words in size. 

After the two files have been properly aligned, a record from each 
file is read using free format reads. The two records are then merged 
by performing a bit-by-bit OR on every common bit. This bit-by-bit OR 
can be described as the logical summing of two bits such that if one or 
both of the two bits have the value 1, then the resulting sum equals 1. 
This summing has the desirable feature of allowing the map outlines, 
whose integer values are 255, to mask out any SASS data jointly located 
in the same spot. This maintains the integrity of the map with only a 
minimal loss to the data set. After the merging is complete, the 
results are written to a new file, and the next two file records are 
read. The final output file is temporarily stored on disk until it can 
be copied to magnetic tape. 
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Figure 3.5 Flow diagram of routine MERGE. 




4.0 PROGRAM OPERATION 

This section is devoted to discussion of the operation of the four 
display generation routines. It is intended to assist those users who 
are familiar with the NOS control system and the Fortran extended 
version 5 language as used by the CDC computer at Langley Research 
Center. The section contains a brief description of the general job 
cards and deck setups. A discussion of each input card and its effect 
on job flow is also included. Reference is made to the source listings 
which can be found in the Appendix. 

4.1 Execution of FALSCLR 

As mentioned earlier, FALSCLR's primary goal is to generate a false 
color display of only the SASS data. It expects input from the data 
tapes and generates a local file containing the display image. 
Intermediate files are created during execution which contain (1) the 
data subset, (2) the raw grid matrix, (3) the smoothed grid matrix and, 
(4) the map projection matrix. Table 4-1 gives the local internal names 
and description of each of these files. 

TABLE 4-1. FILES COMMON TO ROUTINE FALSCLR 


FILE NAME 

DESCRIPTION 

TAPEl 

Un smoothed grid matrix 

TAPE2 

Smoothed grid matrix 

TAPE 3 

Unchanged Mercator Map matrix 

TAPE 4 

Unchanged Polar Stereographic 


map matrix 

TAPES 

Final false color display 

TAPE9 

SASS data subset 

TAPE55 

SASS magnetic tape 


The program can be run in parts by having the proper files available at 
time of execution. If, for example, the smoothed grid file were created 
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at an earlier time, the program could be executed by starting with the 
mapping routines and extracting the needed data from this grid file. 

This option and those like it are controlled by the user's input. 

4.1.1 The Control deck - Routine FALSCLR must be executed as a 
remote batch job because it requires about 127K octal of core memory. 

The source deck is typically on data cards which are read by the card 
reader and the object code is stored on disk using the file name 
FALSCLR. The control deck to perform this operation is given in Figure 
4.1. 

To run the job, the user needs only access the file FALSCLR and 
begin execution. The general job card deck for a complete start to 
finish program is shown in Figure 4.2. In the figure, the 'SASS tape' 
name is the external NOS label of the magnetic tape. The 'user disk 
file name' is any file name the user wishes to give the final false 
color display. Finally, the FOR and EOF refer to the multipunched end 
of record (789) and end of file (6789) cards, respectively. 

To run the program in part, the job card deck must be modified to 
make available the proper data files (see Table 4-1). If, for example, 
the user wishes to only create the data subset and save it on disk, then 
file TAPES is changed to TAPE9 in the deck above. Later, a display 
could be generated from TAPE9 by changing the LABEL card to a 
GET(TAPE9- 'data subset name') and keeping the SAVE(TAPE8=. . . ) in its 
place. Several other options such as this are available by selecting 
the proper file names. Refer to' the flow diagram in Figure 3.2. 

4.1.2 Data Cards - The number of data cards needed for any one job 
varies depending on what operations are done, there are several 
combinations of cards all of which are mentioned here. All input to the 
program uses list directed format where commas and end of cards are the 
value separators. It is necessary to refer to the comments at the end 
of the program FALSCLR in Appendix B to assist in the description of the 
data card deck. 

The first data card for FALSCLR is a program control card which 
sets the starting and stopping points of program execution. This card 
is always the first card in the data deck. It is referred to in the 
input comments as CARD 1. This card detennines the nine possible modes 
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FALSCLR,T50,CM77000. 

USER... 

CHARGE... 

FTN(B=FALSCLR,0PT=2,A,R) 

SAVE(FALSCLR) 

EXIT. 

-EOR- 

source deck 

-EOF- 


Figure 4.1 Creating the task image of FALSCLR. 


FALSCLR,T500,CM130000. 

USER... 

CHARGE... 

GET(FALSCLR) 

LABEL (TAPE55 , NT ,PO=R , F=I , LB=KL , R , VSN= ' SASS tape ' ) 
FALSCLR. 

RETURN (TAPE55) 

REWIND(TAPE8) 

SAVE (TAPE8=' user disk file name') 

EXIT. 

-EOR- 

data cards 

-EOF- 

Figure 4.2 Control deck for execution of FALSCLR 
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of program operation. These modes are combinations of distinct program 
sections and are given as follows: 

Mode 1 - Data Sorting 

Mode 2 - Data Sorting + Gridding 

Mode 3 - Data Sorting + Gridding + Smoothing 

Mode 4 - Data Sorting + Gridding + Smoothing + Mapping 

Mode 5 - Gridding 

Mode 6 - Gridding + Smoothing 

Mode 7 - Gridding + Smoothing + Mapping 

Mode 8 - Smoothing + Mapping 

Mode 9 - Mapping 

Here, mapping refers to the creation of the map projection and 
converting it to the false color file. Depending on the mode, there are 
seven other possible input cards each of which is described in the 

program listing. Table 4-2 shows which input cards are necessary for 
each mode of operation. 

TABLE 4-2. NEEDED INPUT FOR EACH MODE OF FALSCLR 
MODE NUMBER INPUT CARD NUMBERS 

1,2,3, 4,5,6 

1,2,3, 4,5,6 

1,2,3, 4,5,6 

1.2. 3. 4.5. 6. 7. 8 
1,2 
1,2 

1.2. 7.8 

1.2.7.8 

1.7.8 

There are several points worth noting here about certain data 
cards. First, if the grid coordinates are needed (Card 2), then the 
longitudes must be in degrees east. The SASS tapes use a 0 to 360 east 


1 

2 

3 

4 

5 

6 

7 

8 
9 
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scale which causes discontinuity at the Greenwich meridian. The user, 
however, can input a negative value for the minimum east longitude 
whenever the 0 degree line is crossed. The software checks for this 
condition which simplifies the input somewhat. 

Another point to note is that if data sorting is required, then 
multiple modes of data selection are possible. If multiple modes are 
used, then more data cards are needed besides those shown in Table 4-2. 
This may not be clear in the source listing. Card 3 sets the number of 
modes and is followed by that number of data cards defining each mode. 
The observed NRCS measurements will be the same for each mode. 

Finally, if the mapping is performed, the user should recognize 
that the upper and lower data bounds and the number of colors input 
determine the data incrementing value. This value sets the color 
boundary values which are ultimately written on the display with one 
decimal place accuracy. It is recommended that the user check the 
software beforehand to see if the computed incrementing and color 
boundary values are appropriate (use EQS. 2-19 to 2-25). 

4.2 Generating Annotation 

The execution of ANNOT is straightforward and requires only a 
minimum of user interaction. The software is designed so that a 
complete label and map legend are developed from just two input cards. 
The job must be submitted as a remote batch job because it accesses the 
routine RASTGEN which requires up to 177k octal of core memory. The 
final display is usually stored on the user's disk as was done in 
FALSCLR. 

The source deck for ANNOT is typically punched on data cards. It 
is initially compiled and the object module is stored on disk. The 
control cards to perform this operation is shown in Figure 4.3. 

To run the program, the user must access the file ANNOT, the 
graphics package LRCGOSF, and a special library file FTNMLIB before 
beginning execution. Following execution of ANNOT, the user must again 
access the library FTNMLIB, plus the user created routine RASTGEN and a 
user library called SASLIB. These files are used to creat the display 
format image file. The job card deck to generate the display is given 
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ANN0T,T50,CM77000. 

USER... 

CHARGE... 

FTN(B=ANNOT,A,R) 

SAVE(ANNOT) 

EXIT. 

-EOR- 

source deck 

-EOF- 


Figure 4.3 Creating the task image of ANNOT. 


ANN0T,T300,CM177000. 

USER... 

CHARGE... 

GET (ANNOT) 

ATTACH(LRCGOSF/UN=LIBRARY) 

ATTACH ( FTNML I B/UN=L I BRAR Y ) 

ANNOT. 

GET(RASTGEN/UN = ) 

GET(SASLIB/UN=__ ) 

FTN(I=RASTGEN,A,L=0) 

ATTACH ( FTNML I B/UN=L I BRAR Y ) 
FILE(TAPE15,RT=S,BT=C) 
LDSET(FILES=TAPE15,LIB=SASLIB/FTNMLIB) 
LGO. 

REWIND(TAPE15) 

SAVE(TAPE1 5=' user disk file name') 
EXIT. 

-EOR- 

data cards 

-EOF- 


Figure 4.4 Control deck for execution of ANNOT. 
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in Figure 4.4. The 'user disk file name' is again any file name the 
user wishes to call the annotation file. The routines RASTGEN and 
SASLIB are located on another users disk area. Their user numbers have 
been omitted by request, but may be obtained by contacting Frances 
Meissner at the Langley computing facility. 

Upon execution of ANNOT, the user must enter two data cards. The 
first input is a ten character map code used to catalog the final 
display. The second card sets up the color scale annotation and the map 
legend. Refer to the program listing in the Appendix for a better 
description of these two cards. 

The display generated by ANNOT is self-contained, meaning that it 
could byi itself be projected onto the digitizing equipment. The display 
file is somewhat large for an indirect access disk file; however, 
because a particular file could be used repeatedly for displays having 
the same annotation, only a few selections of annotation need be stored 
at any one time. If storage does become a problem, it is recommended 
that the files be stored on magnetic tape and later copies to disk when 
needed. 


4.3 Execution of MAPGEN 

The MAPGEN software has been standardized so that only one input 
card is necessary to generate a unicolored projection of either the 
Mercator or Polar Sterographic map with superimposed land/water 
boundaries. The map projection is scaled before output so that its 
largest dimension is set to a preset pixel size. Because the program 
accesses several large routines, it must be run as a remote batch job. 
Like FALSCLR and ANNOT, the display file that is generated is saved on 
the user's disk. 

The source deck for MAPGEN is usually punched on data cards. For 
convenience, it is read, then compiled, and the object code is saved on 
disk. The job deck to do this is given in Figure 4.5. 


MAPGEN,T50,CM177000. 

USER... 

CHARGE... 
FTN(B=MAPGEN,A,R) 
SAVE (MAPGEN) 

EXIT. 

-EOR- 

source deck 

-EOF- 


Figure 4.5 Creating the task image of MAPGEN. 
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When executing MAPGEN the user must make available the NCAR map 
plotting package and the raster conversion routine RASTGEN. the 
complete job deck to create either of the two map projections is given 
in figure 4.6. Here again, the 'user disk file name' is any name the 
user wishes to give the map display file. The omitted user numbers can 
be obtained by contacting Patricia A. Winters at the Langley computing 
center. It should be noted that TAPE15 (the display file) is copied in 
part to a file called TAPE66 which is then saved on disk. In this 
special usage of RASTGEN, a null file gets appended to TAPE15. This 
null file is undesirable and is eliminated by the copy statement. 

When MAPGEN is executed, one input card is entered which controls 
all the necessary map functions. This card is explained in the comments 
of the MAPGEN listing. Briefly, this card defines the map type and 
geographic region of interest. From this information, a rectangular map 
image is created whose largest side is scaled to 400 pixels. The 
rectangular dimension of the map are output to the user so that they may 
be used as input to FALSCLR to generate the properly sized data map 
projection. It should again be mentioned that the resolution of the map 
image is at present only one degree which must be taken into 
consideration for any map generated by this routine. 

4. 4 Merging the Fil es 

Routine MERGE is a general purpose display file merging routine. 

It expects at least one of the two files it merges to have 512 records. 
The source file is coded on data cards and initially it is compiled and 
the object code is stored on the user's disk. The job deck to do this 
is shown in Figure 4.7. 

Because routine MERGE needs less than 77K octal core memory to run, 
the user has the option of running the program from the real time 
terminals. The option is recommended because it will simplify the 
manipulation of the many files which the user may wish to merge. Also, 
since all working files remain local to the user until he logs off, 
errors can quickly and easily be corrected at one session. The general 
real time control deck is given in Figure 4.8. In the deck, TAPE15 is 
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MAPGEN,T500,CM177000. 

USER... 

CHARGE... 

GET(CLIB,UN= ) 

GET(GLIB,UN= ) 

ATTACH(LRCGOSF/UN=LIBRARY) 

ATTACH ( FTNML I B/UN=L I BRAR Y ) 

ATTACH (TAPE3=MAPDATA/UN= _) 

GET(MAPGEN) 

LDSET(LIB=CLIB/CLIB/FTNMLIB/LRCGOSF,PRESETA=NGINF) 

MAPGEN. 

GET(RASTGEN/UN= ) 

GET(SASLIB/UN= ) 

FTN(I=RASTGEN,A,R,L=0) 

ATTACH (FTNML I B/UN=L I BRAR Y) 

FILE(TAPE15,RT=S,BT=C) 

LDSET(FILES=TAPE15,LIB=SASLIB/FTNMLIB) 

LGO. 

REWIND(TAPE15) 

CO PY ( I =TAPE 15 , 0=TAPE66 , TC=EOF , N=1 ) 

REWIND(TAPE66) 

SAVE (TAPE66=‘ user disk file name') 

EXIT. 

-EOR- 

data cards 

-EOF- 

Figure 4.6 Control deck for execution of MAPGEN. 
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MERGE, T50.CM77000. 

USER... 

CHARGE... 

FTN(B=MERGE,A,R,L=0) 

SAVE (MERGE) 

EXIT. 

-EOR- 

source deck 

-EOF- 

Figure 4.7 Creating the task image of MERGE. 


-log in- 
BATCH, 77000 
GET(MERGE) 

GET(TAPE15='#1 user disk file name') 

GET(TAPE8='#2 user disk file name') 

MAP(OFF) 

MERGE. 

enter data 

SAVE(TAPE44 ='new user disk file name') 
REWIND(MERGE) 

-continue as above or log off- 

Figure 4.8 Terminal Control Commands for Execution of MERGE. 
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the display file which must have 512 records. TAPE44 is the merged 
version of TAPES and TAPE15, and can be saved on disk or merged further 
with another display file. Refer to the flow diagram in Figure 3.1. 

The only input required by MERGE upon execution is the dimension of 
the second display file. If this file has 512 records like the first 
file, a simple one-to-one bit merge is performed on the two files. If 
smaller, offsets are used to properly center the two files before 
merging. In either case, the final file has 512 records. Refer to the 
listing for MERGE in the Appendix. 


5.0 THE DIGITAL DISPLAY EQUIPMENT 


This section is devoted to a brief discussion of the display equipment 
and the steps necessary to operate it. The user who wishes to become adept 
with the system might find this section a useful beginning. Later, only a 
few hands-on sessions are necessary to become proficient. The equipment is 
located at the LRC Computing facility and is available to all contractors. 
Kept with the equipment is a brief manual containing a summary of all the 
steps and commands used to operate the system. 

5.1 Description of the System 

The heart of the digital system is the Intel 80/20 microprocessor. It 
has been programmed to control the various display and tape drive opera- 
tions. The microprocessor has been periodically reprogrammed whenever new 
functions are needed. At present, it has no other purpose than to control 
display generation. 

A CRT is used to communicate with the 80/20. Various three character 
commands can be entered which control tape drive operations, color graphic 
presentation and photographic equipment operations. The CRT is also used to 
display system parameters such as the color scales, the image scales and 
test pattern values. 

The display data is presented to the system by magnetic tape using 
either 800 or 1600 BPI. The tape drive is controlled by a Kennedy Model 
9218 format control unit. Upon command, the Kennedy control unit directs 
the drive to read/write files, rewind or skip several files, and consecu- 
tively read/rewind/read ... a file when that file is being photographed. 

The display image is visually displayed on a 19" Aydin color television. 
This television has an optical resolution of approximately 512 by 460 pixels 
where each pixel corresponds to an eight bit color valve in the tape file. 
The projected image can be shifted around or out of the viewing screen so 
that all parts of an image 1 arger than 512 x 460 pixels can be seen. 

The display image is photographed by a complex instrument called a 
Dicomed Image Recorder. At the option of the user, a photograph can be 
taken of the image on either Polaroid Instant color film or 4 x 5 inch color 
negatives. The recorder performs this task by exposing the film to a 
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miniature CRT which displays the color image one pixel at a time. The CRT 
repeatedly projects the image through a blue, green or red filter to achieve 
the final color photograph. The recorder has three pixel resolutions given 
as low - 512 X 512, medium - 1024 x 1024, and high - 2048 x 2048. For SASS 
display purposes, only low resolution is needed. 

5.2 System Operation 

The digital equipment is usually powered up and ready for operation the 
morning of every vrorking day. If an asterisk appears on the CRT, then the 
system is loaded and ready for use. The first step for the user is to mount 
and load his nine track tape and set the drive to the proper BPI. Once the 
tape drive is online the system is ready to receive data and display the 
image. 

There are several checks the user should make before beginning the tape 
operations. The user should enter the three character commands which 
initializes (1) the color scale the standard test pattern and (2) the image 
scale to the standard size and magnification. The standard color scale is a 
band of colors ranging from blue to green to red to .white in 63 discrete 
steps. This scale can be checked by displaying the system generated test 
pattern. The image size is initialized to 512 x 512 pixels with a one-to- 
one magnification. 

Once the system is initialized, the user should position his tape at the 
desired file by entering the proper tape drive commands. The file size in 
records and eight bit pixels should be known beforehand or can be determined 
by issuing a file check command. Once this file size is known, the user 
must reset the standard image size to a value less than or equal to this 
size. This is needed since the tape drive has no file control to prevent it 
from overshooting a particular file. The data can now be read and displayed 
on the color television. It is noted that the SASS display file uses only 
the lower six bits for color representation. 

The colors used in the SASS display are at the lower end of the 
standard color scale or mainly the blues. This scale will be changed at the 
CRT by first commanding the presentation of the color scale. The user then 
must key in the desired integer to color values. These color values are 
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combinations of percentages of blue, green and red which result in any of 
63^ possible color hues. Some trial and error is needed to obtain an 
acceptable color scale. 

If desired the user can photograph his display by using the imaging 
recorder. The user must first mouht either the Polciroid or negative film 
holders plus film onto the top of the recorder. Next, the recorder must be 
set to the proper resolution and functional color scale. Finally, the user 
must position his tape at the beginning of desired file and then issue the 
proper three character film making commands. The tape will slowly advance 
thru the file three times, where each passage extracts only that data 
containing blue, green or red, respectively. After recording the image, the 
film is removed and developed. The tape file can now be rewound to a new 
file or unloaded if the job is complete. 
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6.0 FUTURE DEVELOPMENTS 

This section is devoted to the discussion of future developments in the 
false color display generation techniques. This report has detailed the 
software as it stands to date; however, it is fully expected that the 
program shall be updated as future needs arise. Several possible 
modifications have already been suggested from the work done, and they will 
be mentioned here. Also several features will be discussed which were 
unnecessary or too complex to be included in the present software. 

6.1 Data Reduction 

As mentioned in section 3.0, the SASS data tapes used by this program 
were simpl if ied versions of the future tapes to be sent from JPL. These 
tapes contain only the basic mode of operation parameters and the 15 NRCS 
measurements. The future SASS tapes will contain a high density of values 
such as (1) multiple cell coordinates, (2) wind speed calculations, (3) 
antenna look angles, and any of a multitude of system identification 
parameters [ref.l]. The software will clearly need to be updated to handle 
this mass of information. 

The primary goal in the reduction of these new tapes will be to create a 
data subset which has only that information necessary to complete the 
desired false color display. This reduction operation will be handled in 
the preprocessing mode of routine FALSCLR and might extract new SASS 
parameters such as antenna incidence angles, wind speeds, and the cells' 
three geographic coordinates. The data subset thus created will then be 
available for any of the various data conditioning functions and ultimately 
presented to the postprocessing routines in the form of a measured SASS 
value and its geographic coordinates. 

It should be noted that for false color displays exhibiting multiple 
satellite revolutions, some sort of data averaging will have occurred in the 
gridding routines in those areas where the revolutions cross each other. 

In order to control this averaging so that the results are meaningful, some 
sort of normal ization must be performed on the data. To date, however, no 
formal method of normalizing SASS NCRS data for off-angle measurements has 
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been established. Since apparently no one particular equation will normal- 
ize the data properly over all measured surfaces, an alternate method has 
been suggested which would apply only to the area under observation. This 
method normalizes the data by first finding the means of the measurements 
from each cell over the area of interest and finally subtracting these means 
from the appropriate NCRS values. This method is being tested and will be 
considered as a future preprocessing operation. 

6.2 Display Generation Changes 

There are several changes that will likely be made in the display 
generation software. Two of the changes arise from the impact of the new 
SASS data tape. The others are changes that have been suggested from the 
results already obtained. This new round of improvements is intended to 
complete the requirement of standardization of the false color display 
techniques for general use. 

The new SASS tapes offer new types and forms of data not applicable to the 
present software. The most significant of these is the three geographic 
coordinates given for each cell. These three coordinates correspond to the 
two cell corners and the cell's center. The cell's single NRCS measurement 
will be assigned to each of these coordinates. The gridding routine will 
need to be adapted to handle these three values such that the two cell 
corner values will contribute less significantly to the final grid matrix in 
relationship to the center value. The grid weighting coefficients for this 
operation will need to be determined before this change can be implemented. 

Several other data types, such as wind speeds and NCRS standard devia- 
tions, might also be presented in a false color display. A few simple 
changes will be needed in the display annotation and scale magnitudes so 
that the program is flexible enough to allow any of these data types to be 
displayed with a minimum of user input. The changes will be implemented as 
the need arises. 

There are two display changes suggested by the results obtained so far. 
First, the map code generated during the annotation is not flexible enough 
to allow repeated use of the annotation files with different data sets. 
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Therefore, a change will be made such that the map code will be 
generated during the final file merging. This code would be a four 
digit number which would uniquely define each display. 

The second display change would be to print the latitude and 
longitude values on the map image, something the present software does 
not do. This improvement will probably have to be accomplished through 
the use of some of the subroutines available in the SUPMAP Graphics 
package. The final coordinate format will be determined when it has 
been ascertained which graphics routines can be used. 

Finally, it has been suggested that other map projections be 
incorporated into the software. The most desirable projections, besides 
those already in use, would be the Lambert Conformal projection. This 
map type is ideal for the midlatitude regions [ref. 3]. The SUPMAP 
package contains this projection as one of its nine map types. To make 
use of this option, a grid-to-Lambert projection routine will be 
required in the postprocessing mode of FALSCLR. The software for this 
routine is in the beginning stages and will be tested and implemented 
upon completion. 
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^PPCG.AM "ALSCLR(INPUT.^U;PUT^TAgE3=^IN^^^^ 

THIS PROGRAM WAS CREATED ON APRIL 12. 1979 BY CLAYTON B. JACKSON WHILE 
**♦ UNDER THE EMPLOYMENT OF THE RESEARCH TRIANGLE INSTITUTE OF N.C.. IT til 

*** rpNTOArT'"wf§i CONTRACT REQUIREMENTS UNDER NASA *** 

NASI-153S8. THE SOFTWARE HAS BEEN REVISED AND TESTED AS OF ♦♦♦ 
JUNS 19* 4(4(1^ 

PPOGRAM IS ^DESIGNED TO GENERATE FALSE COLOR DISPLAYS OF SEASAT-A 

WITHIN SELECTABLE GEOGRAPHIC REGIONS. IT IS INTENDED 
CCNTPOL CATA*S CYBER 170 SERIES. CYBER 70 SERIES. SERIES 
6000 AND SERIES 7000 COMPUTERS USING THE FTN5 LANGUAGE. THE PROGRAM 
UNIT PERFORMS THE FOLLOWING FUNCTIONS: PHUGRAM 

*“ SORTED BY GEOGRAPIC LOCATION AND VARIOUS SASS 
PARAMETERS. AND THE RESULTING DATA IS WRITTEN TO TAPE9. REFER 
TO POUTINE DATAIN« 

2) . above IS READ AND THE APPROPRIATE DATA IS SELECTED FOR GRIDOING. 

the FINAL GRID MATRIX IS WRITTEN TO TAPEl. REFER TO ROUTINE SSGRID. 

3) . THE GRIDDED DATA (TAPEl) IS SMOOTHED BY INTERPOLATION FOR MISSING OR 

• THIS SMOOTHED GRIDDED DATA IS WRITTEN TO TAPE2. REFER TO 
hCUTINE smthcd* 

A). THE GRIDDED DATA (TAPE2) IS PROJECTED ONTO A MERCATOR ANO/OR POLAR 
ST FREOGR APH I C MAP, THE FINAL MAP DATA VALUES ARE INTEGERS 
PFoRESENTING FALSE COLOR SCALING OF THE SASS DATA. THIS MAP DATA IS 
WRITTEN TO tapes, REFER TO ROUTINES MERCAT, POLAR, COLORS AND PACKEM. 


^ FUNTTONS NEED NOT BE RUN CONSECUTIVELY. ANY OF THE FIRST 3 
STEPS COULD HAVE BEEN PUN PREVIOUS AND THE RESULTS STORED ON DISK OR 
MAG tape, various INPUT CARDS ALLOW YOU TO ENTER AND LEAVE THE MAIN 
PROGRAM AT SEVERAL LOCATIONS. 


***** PLEASE REFER TO END OF MAIN PROGRAM FOR COMPLETE DESCRIPTION OF ALL 
THE NECESSARY DATA CAROS ***** 


D I MENS ION GRID<40D.20),SUM(A00,20). INDEX! 40l),P{A00),SS<400) . 
* I FILL (400 ) 

CCMMCN/XMAP/ULT.DLT.ULG.DLG 
MAXS7=400 
PRINT 3999 , 1 

C 

C ♦* PROGRAM CONTROL INPUT 
C 

RFAC(5,*) IGPIO, INTERP, IDATA 
IF(IGRID) 135,4E,45 
45 PRINT 3999 . 2 
C 

C ** GRID COORDINATES ENTERED 
C 


REAC(5,*) ULT.DLT.ULG.DLG, GROSZ 
IF< DLT.LT.-35. ) DLT=-65.0 
IF(ULT ,GT .85. ) ULT=a5. 

C 

C ♦* SET GRID MATRIX DIMENSIONS 
C 

M=( ULT-DLT) * U 1 . I 1 2/GRD3Z*! . 

I F ( M ,GT .MAXSZ ) GO TO 399 
IFICLG.LT.ULG) GO TO 21 
DL G = DLG-360. 

21 CCNTINUE 

N=(ULG-DLG)*1 1 I . 1 12/C,RCSZ+1 . 

IF ( N .GT .MAXSZ ) GO TO 399 
IF( IDATA) 155.175.155 
C 

C ♦* CALL DATA SELECTION ROUTINE 
c 

175 CALL DATAINI IDATA, 55) 

1F( IDATA. EC. 0) GO TO 699 
ISE REWIND 9 

IF( INTERP. LT.O) GO TO 599 
IF( IGRID.GT.4) GO TO 255 
^ , WRITE! 1 , 302) N.M, ULT.DLT.ULG.DLG 

C ♦* PERFORM GRIDOING 
C ■ 

CALL SSGPIC(N.M.ICATA.GRI0.SUM.20) 
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REWIND 1 

IF( INTERP.lt .2 ) GO TO SS 
255 CONTINUE ^ ^ 

WRITE<2.302> N.M.ULT.DLT.ULG.DLG 
REACH .302) 

*♦ SMOQTR GRID MATRIX 

CALL SMTHGDI INTERR. GRIO.N.M .20.SUM) 
REWIND 2 

IFI IGPID.EQ.O) GO TO 499 
GO TO 35 
READ 11.302) 


55 

♦ 4 


WRITE FILE TAPE2 


WRITE (2. 302) N , M . UL T , DLT . ULG. OLG 
DO 36 1=1 .M 

READ(I.401) (P(J).J=l.N) 

36 WRITE(2,40l) {P(J).J=1.N) 

ENDFILE 2 
ENDFILE 2 

rewind 2 

IFI IGRID.EQ. 0) GO TO 499 
WRITEI6. 4003 ) 

GO TO 35 

135 PEAD(2.302) N. M . ULT .DLT .ULG .OLG 
REWIND 2 
'print 3999 . 7 


♦♦ MAP TYPE AND DIMENSIONS ENTERED 

35 PEACI5.*) N1.N2.N4. ZLAT.MAPTYP 
IFIEOFISn 40.25.40 
2S IREC = N2'*’ I 

I F( M .GT .MAXSZ.OR .N2.GT .MAXSZ ) GO TO 199 
IFI MAPTYP.GT. I ) GO TO 31 

♦ 4 PEPT'CRiM MECATOR MAP PROJECTTON 

CALL MERCATIN.M.SS.N1.N2.1NCEX.IFIUL.P.IREC) 
ENDFILE 3 . „ , 

WRITEI6.4001 ) Nl .N2 .ULT .DLT .ULG.OLG 
GO TO 40 

31 IFI N4.EG.O.OR.N4.GT .20) N4=20 

B 1=N2 

32 B2=B1/N4 

1FIAINTIB2) .EG.62) GO TO 33 
N4=N4- I 
GC TO 32 


44 

33 


PERFORM POLAR STEREO uR APH I C PROJECTION 


CALL POLARIN.M. ZL AT. GRID. Nl . INDEX .P.IFILL.N4.IFEC) 

ENDFILE 3 

WR ITEI6.4002) Nl . N2 .9 0 . . ZLAT . 360 • . 0 . 

40 ENDFILE 8 
REWIND 8 
REWIND 2 
REWIND 9 
STOP 

399 WR I TEI 6 . 3000 ) ULT .DLT .ULG , OLG . GROSZ 
STOP 

199 WRITEI6. 3001) N1.N2.. MAXSZ 

STCP 

499 W RIT E I 6 . 3 002 ) N . M . ULT . DLT .ULG . OLG. GROSZ 
REWIND 9 

I FI 1NTEPP.lt .2 ) GO 501 

WRITEI6.3003) INTERP 
STO = 

SOI WRITEI6.3004 ) 

STOP 

599 WRIT5I6. 5000 ) 

STOP 

699 WPITEI6.500 1 ) 

STOP 

302 FCRMATI 21 3. 4F1 0. 4 ) 

3000 F CPM AT 1 1 HO . 5 0HTFE AREA OF INTEREST IS TOO THE GRID SIZE 

Wl'^H IS TOO SMALL— —. X 1 X . 35HCHANGE THE FOLLOWING TO GET NO MORE. 
422H THAN A AOO X 400 GR ID . / / 1 X . 5HNL AT= . F 7 , 2 » S X . 5H SL AT- .F7 . 2 . 5X . 

47HELCNGl=,F8.2.5X,7hELONG2=.F3.2.5X .lOHGRID SI ZE=.F6. I .4H KM. ) 

3001 FCRMATI IHO. IOHEITHER N 1= . I 4 . 2X , 6HQR N2= . I 4 .2X, I 2H A RE GREATER , 
45H THAN.T4.2X.39HREENTER NI.N2-THE FALSE COLOR GRID SIZE) 

■>nnp FORMAT! 26H0THE CCMPLTEO SASS GRID IS.I4.10H PIXELS 8Y.I4. 

45 1 H PIXELS AND 19 LOCATED WITHIN THE FOLLOWING REG I CN I . // 1 X . I OX . 

^liNLAT=.F|.2,5X.5HSLAT=.F7.2.5X.7HFLONGl=.F8.2.5X.7HFUONG2=.F8.2 

4//1X. IOX.36HANO the GRID BOXES ARE APPR OX I MATEL V ,F7 . 2 . 1 X . 
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nonnnnnnnnnonnnnnnnnnnnnnnnnnnoonnnonnonnnnrinonnnnrinonnnnnnnnonririrtnnnoo 


♦ <5HKM SQUARE) 

3003 FOPMAT(51HOTHE 1 0 WAS SMOOTHED BY INTERPOLATING FOR MISSING, 

♦ 30H DATA USING AN INTEPP VALUE OF , 1 2, 2X . 1 7H ( REFER TO SMTHGD)) 

300A FORMAT ( 40HONO SMOOTHING WAS aERFCRMED CN GRID DATA) 

3P99 FORMAT ( I 6H0ENTEP CARD NO. ,!l) 

4001 FORMATISIHOA MERCATOR MAP PROJECTICN OF SASS DATA WAS CREATED, 
*18H WITH GRID SIZE OF.I4.IOH PIXELS BY,I4,7H PI XELS ,/ 1 X ,6HAND IS. 

♦ 33H LCCATED IN THE FOLLOW ING REG ION 1 . // 1 X , I OX , 5HNL A T= , F 7 . 2 , 5X , 
45HSLAT=.F7. 2 , 5X , 7HEL0NG1= ,F8 . 2 , 5X , 7HELC NG2= , F8 . 2 ) 

4002 FORMAT(54HOA POLAR ST EREOGR APH I C MAP PROJECTION OF SASS DATA WAS, 
*26H CREATED WITH GRID SIZE CF,I4,10H PIXELS 8Y,I4,7H PIXELS, /IX, 

♦ 39HAN0 IS located IN THE FOLLOWING REGI ON : , // I X . 5HNLAT= ,F5. 1 . 5X , 
*5HSLAT=,F7.2.5X.7HEL0NG1=.F6.I , 5X . 7HEL0 NG2= . F 4, I ) 

4003 FORM ATI IH 1 ) 

5000 FORMAT (47H0CNLY ROUTINE OATAIN WAS PERFORMED CN SASS DATA) 

5001 FORMAT(39HONO DATA WAS FOUND IN AREA OF INTEREST-, 

♦17H CHECK DATA CARDS) 


THE FOLLOWING IS A COMPLETE SUMMARY OF EACH INPUT CARD 

CARE MUST BE TAKEN WHEN WRITTTNG THE INPUT CARDS SINCE THE PROGRAM 
DOES NOT CHECK FOR IMPROPER INPUTS. 

** ALL INPUT IS LIST DIRECTED WITH COMMAS USED AS FIELD SEPERATORS ♦♦ 
THE CARDS ARE AS FOLLOWS: 


CARO I. ENTER PARAMETERS IGRIO. INTERP AND ICATA WHERE: 

( <0 PPQCFDE DIRECTLY TO MAPPING 

IGRID= ( =0 PERFORM ALL BUT MAPPING 

( >0 <5 PERFORM EVERTHING 

( >4 PERFORM SMOOTHING AND MAPPING ONLY 

( <0 PERFORM CNLY DATA SORTING AND SELECTING 
INTEPP=( =t PERFORM EVERYTHING EXCEPT SMOOTHING 

( >l PERFORM EVERYTHING- INTERP SETS AMOUNT OF 
SMOOTHING. REFER TO ROUTINE SMTHGD. 

( <0 DATA SUBSET DONE EARLIFP- UNKNOWN AMOUNT OF DATA 
IOATA= < =0 PERFORM DATA SORTING AND CONTINUE 

( >0 DATA SUBSET ALREADY OCNE- IDATA=NO. DATA POINTS 


CARD 2. ENTER GRID COORDINATES - NORTH AND SOUTH LATITUDES AND MOST EAST 
AND LEAST EAST LONGITUDES- ALSO ENTER GRID BOX SIZE IN KM. 

THIS CARO OMITTED IF I GF I D IS NEGATIVE *** 


4 . 4 * 444 * OMIT CARDS WITH ASTERISKS IF IDATA IS NONZERO CR I GR I D < 0 ****** 

444 CARD 3. ENTER NMOOES- THE NUMBER OF DATA SORTING CONTROL TRIOS. 

*4 IF NMOOES FGUAL ZEPOINC SORTING) CMIT CAROS 4 TO 3HNMODES *** 


♦** CARD 4. ENTER A SASS MODE, BEAM NUMBER AND POLARIZATION ON EACH OF 

4*4 CARO 5. THE NMOOES INPUT CARDS. ENTERING ZERO FOR ANY PARAMETER 

*44 , . CAUSES NO SORTING FOR THAT PARAMETER FOR THAT PARTICULAR 

4*4 . . INPUT CARO. 

♦4* CARD 3HNMCDF5 


4*4 ASSUMING NMOOES=l (FOR EXAMPLE) THE INPUT CARDS CONTINUE AS: +** 

*44 CARD S. ENTER CELL NUMBERS TO BE USED IN SORTING-44 A BLANK (ALL 0*S) 
CARD IMPLIES ALL CELLS ARE USED. 


444 CARO 6 . ENTER CCOPDINAtes cF DATA SUBSET- NORTH AND SOUTH LATITUDES AND 
MOST EAST AND LEAST EAST L C NG 1 T JOE S . 4 4 NEGATIVE LONGITUDES MAY 
BE USED EVEN THOUGH DATA TAPES 00 NOT. 44 


44*444 OMIT ALL FURTHER INPUT CAPOS IF IGR 1 0= 0 OR INTERP < 0 444444 

CAPO 7. ENTER MAP SIZE IN PIXELS- WIDTH X HEIGHT- ENTER MAP WORKING 

HEIGHT <=20 PIXELS - EN'ER ZLAT, THE SOUTHERN MOST LATITUDE FOR 
POLAR STEFFCGF APHl C MAPS - AND ENTER M4PTYPE WHERE: 

MAPTYP"=( <=1 MERCATOR MAP PROJECTION 

( >l POLAR STEREOGRAPhIC MAP PRCJECTION 


CARD 3. ENTER FALSE COLOR PARAMETERS ISCALE. ICOLCR, IBQUNO. SIGHt, 
AND S ICLO where; 
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nonooaoooonnnno 


I sc AC e= 


(I LINEAR COLOR TO DATA PELATIONSHIP 
(2 LOGARITHMIC COLOR TO DATA RELATIONSHIP 
(3 EXPONENTIAL COLOR TO DATA RELATIONSHIP 


tCOL0F= NUMBER OF FALSE COLORS 01 SPLAYED <= 20 

ie0UND= (I DATA ABOVE AND BELOW LIMITS IS NOT COLORED 

(2 DATA ABOVE AND BELOW LIMITS BOTH HAVE A COLOR 


SIGHI = MAXIMUM SASS DATA VALUE TO BE SCALED 
STGLO = MINIMUM SASS DATA VALUE TO OE SCALED 

END 
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SII8!»0rj?INE DATAIV (IDATA.rUNIT) 

THIS ROUTINE ESTABLISHES A SASS DA” A SHH^ET AMO rnn-pflTs I'T TO PtlP TAPE‘5. 
THE oppRRjM TJS5S TUP PEESTftT fci’EA”' AMD SOt^S "’H’^ S^SS DATA PI US'^H 

SELFC^ED MaiTIPLES OF ECDES- BEA’* MIJMBFFS AMD DO L A tI E ATIO NS FOB ANY OP ALL 
OF '”HF ANTENNA’S 15 CELLS. THE G”''C'’A?HIG bEOION OM Tfp SB3SFT is spi 
WITHIN THIS PCUTINE AND DOES NOT MECFSSA'JI.Y rro t xir) to Tpp GPIDDED 
APFA. EEFEP ""O NAIM PECGFAN FOE D 3 SC « T dtt r; N op ”NPPt CARDS. 


DIMENSION MODES (1 5 ) , I B EA MS ( 1 5) , I P^I S f 1 5 ) , TC ELLS ( 1 S) 

PEAL NLAT ^ i i 

tdati=o 

ASSIGN 1 TO IGO 
ASSIGN 9 TO JGO 
ASSIGN 12 TO MGO 
PRINT 3999 , 3 
READ(5,*) NMODES 
IS(nnODES) 20,20,35 

♦ » DATA SFLECTORS MO OE , BEA M , POLA P 17, A” T ON AND CFT.L NOS . INPIJT'”ED 


35 

130 

20 


*♦ 


DO 110 I = 1,N!10DFS 

PRINT 3999 4 


P.FAD(5,*) MODES (I) ,IBFAMS (T) , 
PPIN"’ 3999 ,5 
FFADf5,*) (ICELLS (J) , J= 1, 15) 
IF(ICELLSJI) .EO.O) AsSIGN 13 
PRINT 399^ ' f, • 


I POLS (T) 


TO MGO 


** DATA SHBSET COORDINATES INPUTTED 

PFAD(5,*> NLAT,SLAT,ELGHI,ELGLO 
XPLG=ELGLC 

♦♦ CHECF FOR NEGATIVE LONGITUDES 


IF (ELGLC.lt. 0.] 


IF(ELGLO.LT.O. 
IF (NLAT.GT. 


XDLG=360. aFLGLO 
ASSIGN 2 TO IGO 

. , . . ASSIGN 10 TO JGD 

IFfSLAT.GE.HLATl ASSIGN 1C TO JGd 
I F (FLGLO.lt. ELGHI) GO TO 56 
if(elglo.eo.o.) go to 156 
XDI,G=PLGLO 
ELGLO=ELGLO~36C. 

ASSIGN 2 TO IGO 
GO TO 56 

156 ASSIGN 10 TO JGO 
56 CONTINUE 

** read SASS TAPS 

150 READdUNIT, 200) MODE , IBEA M , I POL , A ZV SA T 
I? (EOF (ItJNIT) ) 30,132,30 
132 IMFASpO 

IF (NMODES. EO.O) GO TO 170 


PEPFODM DATA SELECTION 
DC 14(J 1 = 1, NMODES 

-j 

GO TO 140 
GO TO 140 


I i — I « IN 1 0 ; IJ 12 o 

IFfMCDES(I)) 4,4,3 
IF hODF.NE. MOOFS fl) ) GC 
IF (IPEAMS (I) ) 6,6,5 

IF (IDEAM.NE. tBFAMS (I) ) 


IF(IPOLS(I)) 8,8,7 
IFjlPOL. NE. TPOLS (I) ) GO TO 140 
IMEAS=1 
140 CONTINUE 

IF(IMEAS) 170,33,170 
1 70 DO 120 1=1 ,15 

READdUNIT, 201) TO ELL, SIGZ , XLA t’ , XLO N’G 
IF (EOF (lUNIT) ) 30,133,30 

133 GC TO "GO, (12, 13) 

12 IMEAS=0 

DC 180 J = 1, 15 

IF (ICELLS (J) .EO.O) GO TO 15 
IF (ICELL.NE. ICELLS (J) ) GO tg IPO 
r."EAS = 1 
180 CONTINUE 

IFdHFAS) 13,120,13 


15 

1 3 

2 
1 
9 


GO '”0 IGO, (1,2) 

IE (XLONG.GE. XDLG) XL0NO,= XL0NG-3 60 . 

GC TO JGO, (9, 10) 

IF (XLAT.GT.NLAT.OR.XLAT.LT.SLAT) GO 120 
IF (XLONG.GT. ELGHI.OP. XLONG. LT. EI.GLO) GG PC 

WRITE FILE TAPE9 CONTAINING DATA SUPPEt 


12 ' 
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Of Id 



WFTTE {6 ,400) IDATA 
RE^fJPN 

101 FCRMAT(I4). 

102 FCP1A'^i3I3» 

103 FCPflAT 1513) 

104 FCPIATl 4F10.0) 

200 ?CR!1ATi 3I3,F10.2) 

201 FCPf1A?(I3,4rg.2) 

3C1 PPPFA? (3F10.4) 

400 '='0 PH, i\ >7 (31 HO** MO. OF OOIDDFO DAT." P''T M"'3= , I fi , PR **) 
3994 FCPHAT(16HCFNTFP CARD NO. ,11) 

HMD 
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c 

c 


c 

c 


120 

c *• 

51 


S MBFOriTIME SSG3 1 D (•! , M , ND , (IP ID, SO v , 2) 

GPTDS THE SASS DATA A IPTD EAT^tk WHEF'’ T’HF fJRID 30YES 

user iNPnrrED size cf yxy pv,. all sass data fpom thf self^^’^d data 

FALL INTO EACH IKDIVIDHAL GPTD B'^X f c; ' A V FP A G ED BX '^NVFRiip 
SPUAPE DISTANCE WEIGHTING AND IS' SIPPED IN THE FINAL GPIDHATPIx!^ 

FEAL GPID(N,H2) ,Sn.*I (N,N2) 

Cf^'1NON/XMAP/ULT,DLT,niG,DLG 
DATA BAD/-8B.8a/ 

NP = ND 

IF(ND.LT.O) ND=5C0C0 

AN = N 

AN2=N2 

NX=0 

ASSIGN 1 TO IGP 
IF (DLG.GF.P. ) GC TO 120 
IDI,G = DLG>360. 

ASSIGN 2 TO IGO 
CON'^imiE 


SET GRID BOX DINENSIONS 

DSEG1= f(JLG-DLG) /»! 

DS^G2= (tILT-DLT) /M 
NX = NX«-N2 

SENT NO 9 

,7HP=N2 

C *♦ IN"’IALIZS GRID MATRIX 


30 


10 

15 


26 


22 

20 


c ** 

70 


DO 30 1=1, N 
DC 30 J=1,N2 

sn« (I,J) =i5. 

GPID (T,J) =0. 
rLAT=0t.T -DSEG2*MX 
UIAT=ULT“DSFG2* (MX-N2) 


PF AD DATA FILS 
DO 10 1=1, HD 


'APE9 


25 


I II i= I , HU 

PFAD(9,100 SIG,XLAT,'fLONr 
IF (EOF (9) >15,25,15 
OOtHiNfJE 


c 

c ♦* 

c 


f F (SIG. LE. BAD. OP. XL AT. FQ. BAD. OP. XLONG, B AD) G''^ 10 

GO TO IGO, (1,2) ' ' ■ 

IF (XLCNG.GT. XDLG) XLONG=XLONG- 3 60 . 

CCNTINDE 

IF(XLAT.GT.ULAT.OR.xLAT.L'".DLA'^) "O 10 
IF (XLCN3.GT.tJLG.0R. XLONG. LT.DLG) GO TO 1-? 

COMPOSE GRID VALDES 

Bl= (XLONG-DLG) /DSFG1 
B2= (DLAT-XLAT) /DSEG2 
I1=MIN1 (AN.B1+1) 

I2=MIN1 (AN2,B2+1) 

P2=(B1-T1 + .5)**2i-{B2-I2+.5)*i‘2 

IF (R2. LT. 1 . E-B)^ R2=1.Sr-8 

3 '’ID (11,12) =GRID (11,12) 1-SIG/P2 

SIJM(fl,I2)isDH (li,l5) ♦1./R2 

CONTINDE 

Nn=I-l 

DC 20 1=1, N 

DO 20 J=1,N2 

I 22,22,26 

GRID (T,J)=GPID (I,jf/SUM(I, J) 

GO *00 20 


C ♦* FLAG MOLL GPID BOXES 


6 1 

100 

300 


GPID (I, J) = 1. E5 
CONTINDE 

IF(MX.GT.M) J0P=N2l-M-!1X 

NPITE GPID MATRIX TO FILE TAPEX 

DO 70 J=1 JOP 

WPITE(1,30C) (GRID (I,,J) ,T=1 , N) 
TF(MX.GE.M) GO TO ^ 1 ' ' ' ' 

GO 51 

ENDFILE 1 
FNDFILF 1 
FORM AT (3F1C. 4) 

FORMAT (10F1 1 .4) 

D P N 
END 
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c 

c 


?tIDPOrjTINE S’fTaoO (IG,3,N,>1, M2,XSnt1) 

THIS POOTINE INTERPOLATES S ASS DAT? VAI.U 
VALUE DOE TO MISSING DR SPARSE DA""?-. A M 
THE 5X5 GRID BOX HATPIX SORROONDTMG TH 
DETERMINED BI THE INVERSE SOOARED DISTAN 
E*inTY GRID BOX. PAEAHETEP IMT’^RR IS 

ITS VALUE DETERNTNES THE AMOUNT OF GFI'a 
5X5 GRID MATRIX WHERE THE TSNTF'’ BOX I 
THF BOXES AP.'’ POSSIBLE VALDES OF INTRRP. 
LESS THAN OR EQUAL TO THE INPUTTED VALUE 
THE SHCOTHING PROCEDURE. 


_S FOR GRID BOXES WHICH HAVE NO 
"TOKTED averaging IS DONE FROM 
’JMPTY GRID BOX. THE WEIGHTS ARE 
OF from ■OitE KNOWN GRID BOX TO THE 

TMp(j"">rn TM 7 he main program and 

3M''OTHING. BELOW IS A SAMPLE OF A 
S "’0 BE SMOOTHED. THE NUMBERS IN 
THEREFORB, ""HE BOXES WITH NUMBERS 
OF INR’EPP ARE THE ONES USED IN 



C 

c 

c 

10 

57 

U5 


U4 

C 

0 ** 

c 


35 


34 


33 


32 


31 


2B 

29 

30 
C 

c ♦♦ 


REAL G(N,N2) ,XSUM (N,N2) 

D»R'A BtAG/1.E5/ 

DATA W1,W2,W3,W4,W5/1. ,0.5,0.25,0.2,0.125/ 

IN0LL=1 

IF(IG.GT.3) INULL=2 

N3=INULL*2+1 

ISP=N-INULL 

JSP=M-INULL 

JST=INULL 

IST=INULL«-1 

READ GRID FILE TAPE1 

DC 10 1=1, N3 

WRITE(2,100) (G(K,I) ,K=1,N) 

J=IST 
JST=JST+1 
DO 44 I=1,IN0LL 
XSOM (1,1) =G (I,. I) 

IXZ=TSP4-I 

XSOM (IXZ, 1)^=G (IXZ, J) 

DO 30 I=IST,ISP 

compute SMOOTHED GRID VALUES FOR MISSING DATA ONLY 

IF (G (I, J) .NE.FLAG) GO TO 2S 

DSUM=0. 

rcT=o. 

GO TO (30,31.32,33,34,35) ,IG 

SIIM=G (t-2,J-2).+G (T-2,0+2) tG(I + 2,J-2) fG(If 2, J + 2) 

ITOTO= (S0M*5. E4)/FLAG 

nSUM = DS0M<- (SUM-ITOTO*FLAG) *W5 

TCT=T0T+ (4. -ITOTO) *W5 

SnM = G(r-1,J-2) +G (1*1, J-2) +G (I-2.J-1) + G fI+2,J-'') 

»+G (I-2,Jf1) +GJ[T + 2,J+1) ♦G(Im1,J«-2) +G (T+1,.1 + 2) 

ITOT0= (S0M + 5.E4)/HaG 
DS0M=DSUM* (SUM-TTOTO*FLAG) *W4 
T0T=T0T* (B. -ITOTO) *W4 

SUM = G (I,J-2)«-G (1-2, J) *G (I4-2,J) ♦Gd, J + 2) 

XTOTO= (SUM+5. E4) /FLAG 
DSUM=nSUMA- (SUM-IT0Tn*FLAG) *W3 
TOT=TOT?- (4.-ITOTO) *«3 

SUM = G (I-ljJ-lf+G fl+1,j-1) to, (1-1 ,J+1)+G (H-1 ,J+1) 
TTOTO= (SUMtS. EU1/FLAG 
DS0M=DSUH+ (SUM-TTOTC*FLAG) *W2 
Tf''T>=TOTt (4.-iroT0) *W2 

S0M=G(T,J-1) to, (T-1 ,J) tG (I+1,.1) +G 

ITOTO= (SUMtS. S4) /FLAG 
DSUM=DSUM+ (SUM-ITOTO*FLAG) *W1 
T'G'r='roT+ <4 . - ) * W 1 

IF(TOT)' 29,28,29 
XSOM (I,1)=G (I, J) 

GC "BO 30 

XSOM (T,1)=DSUM/TCT 
CONTINOE 

WRITE SMOOTHED GRID MATRIX TO FILE R'APF2 
WRITE (2, IOC) (XSOM (II. 1) .11=1 ,N) 
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ub 

5C 

55 


10? 


DO 5D K=2,(J3 ' ■ ■ ■ 

DO, 50 1 = 1 In 

n(t,K-i)=(5(r,K) 

PFAD(1,100) (0 (II,N3) ,II=1,N) 

IF 50Ff1) )' 55,46,^5 
!!4 = N3-1 

DO 50 I=TST,M4 

»PI”|^2,100) (G(TI,I),ir=1,H, 

ENDFILE 2 
FFWIND 2 
PEHIND 1 

BrT'(PM 


FCPMAT(10F11.4) 

END 
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SUBROO?INS HBPCar (N,fI,X,N2,N1,TND''X,I?ILL, Y,IP) 

THIS BPOTINE PEPFOPHS A PROJECTICX OF SASS ROIDDED DATA ONTO A MERCATOR 
MAP HAVING THE SAME GEOGRAPHIC COORDINATES. THE OSER SHPPLIBS THE SIZE OP 
THE RECTANGULAR MERCATOR MAP MATRIX IN PIXELS DP TO HOO X 405 PIXELS. 


DIMENSION INDEX (IR) ,X (N) ,X(N2) , IFILL (N2) 
CCMHON/XMAP/ULTtDLT.HLG.OLG 

DATA A1 ,A2,A3/.7853§6 1635, R. 726 646 2609F-3, 1.745 3292 52E- 2/ 
CALL SECOND (TX) 

HRITP(6.500) 1,TX 
READ(2,201) 

CALL OPENMS (3 , INDEX ,IR ,0) 

(A1tOLT*A2) ) 

Ul + DLT*A2) ) 


AX=ALOGfTAN 
BX=ALOG (TAN 


PROJECT EQOATOR LOCATION 

E0DAT=N1*AX/(AX-BX) 
SCALE=EQOAT / AX 
XLAT= (DLT-DLT 
XINC=1.*N2/N 
JSTR=1 

DO ICO K=1,H 
ISTP=1 


/ H 


** READ GRID PILE TAPE2 
PEAD(2,101) X 

*♦ PROJECT GRID VALDES TO MAP MATRIX 

DO 10 1=1, N 
ISTP=I*XINC+0. 5 
IF (ISTP.LT.ISTR) GO TO 10 
DC 20 J=ISTR,ISTP 
20 Y(J)=X(I) 

ISTP=ISTP+1 
10 CONTINUE 

DLAT=DLT-K*XLAT 
CX = AL0G (TAN (A 1+DLAT*A2) ) 
rSTP=BQUAT-SCALF*CX 
IP (ISTP.GT.NI) ISTP=N1 
IF(ISTP.LT. JSTR) GO TO 100 
DO 60 L=JSTR,ISTP 
CALL HRITMS(3,y,N2,L,0,0) 

CONTINDB 
JSTR=ISTP+1 
CONTINOE 

IF(ISTP.EQ.NI) GO TO 33 
DO 34 I=ISTP.N1 
34 CALL «PITMS(3,y,N2,I,0,0) 

33 CALL SECOND(TX) 

CALL SECOND{TXf 
WRTTE(6,500) 2,TX 

*• INVOKE COLORS ROOTINE 

CALL COLORS (N1 ,N2, X ,IFILL, INDEX , 3 , 1 R) 

CALL CLOSMSny 
CALL SBCOND(TX) 
aRITE(6,500) 3,TX 
101 F0RMAT(10F11.41 
201 FOPMAT(2I4,4F10.4) 

500 FORMAT (9HOSTEP NC.,I3,18H IS ROUTINE MERCAT,10H CCM PLETED, F8. 2 , 
♦36H SECONDS AFTER START OP MAIN PROGRAM) 

RETOPN 

END 


60 

100 
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c 

c 

c 

c 

c 


c 

c; ** 

c 


2 


C *» 

c 


c 

z *« 

c 


10 


3C 

c *♦ 

4C 


m 


21 


1 1 


13 

c 

z ** 

c 

20 


SOBPCtJTINS POLAP (*1,M,ZLA'^,DS,:n ,r''nS'(,o,T5'rLL,VU,r3) 

PHIS POaPINF PEPFOPKS A PPCJFCmH PASS G"inDFD DATA ONTO A POLAP 
STEPECGPAPHIC MAP HHOS GE'^GRAPHIC CEVT?P IS THE POLE AND WHCS SOUTHERN 
M^ST lATITHDE IS DSER INPUTTED. TH'’ MA” PPriJ’^CTTON IS A SO'JAPE MATRIX 
HHEPE DSER SOPPLISS THE SI3E IN P’-XEIS HP TO 4f 0 X 400 PIXELS. 


DIM’" NS ION INDEX (I P 1 . P f N 1 ) , DS {N . 
,DNLT,0L0,DL3 


, N4) , I’^TLL (N 1) 

C^MMCN/XM AP/HPLT,DNLT,OLfi,DLG ' 

DATA A1 , A2, A3/. 7H5 3 Q« 1 6 3S , 8 , 72 6 64 fc 2f C ^ 3 , 1 . 74F 3 29 2 52E- 2/ 

A4=57. 29577961 

FLAG=-8H.88 


SET MAP TO N^FTH OR SOUTH POLE 
PI.AT=ZLAT 

IF (UPLT. LE.C. . AND.DNLT. LT.C.) PI.A' 
ULT=AMAX1 (ABS (UPLT) , ABS (DNLT) ) 
DLT=AMIH1 (ABS (UPLT) ,ABS (DNLT ) 
AH3=N1/2.tC.S 
CALL SECOND (TX) 

H?ITE(6,5C0) 1,TX 
PEAT (2,601) 

N3=N1/2 
DO 2 1=1 ,N1 


'=-c>I,AT 


P (I) =FLAG 
CALL 


OPENMS (4,INDSX,IR,0) 


SET 3PI0 DIMENSIONS 


XLAT= (ULT-DLT) /M 

XlONG= (ULG-DLG) /N 

SCALE = AH3/TAN (A1-PLAT«A2) 

INITIALIZE MAP MATRIX FIIE 

DC 10 1=1, N1 

CALL WRITMS (4, P,N1 ,1,0 ,C) 

CONTINUE 
CALL SECOND (I’X) 

HFITE(6,5D0) 2,TX 
HLAT=ULT-90. 

P1 = SCALS*TAI! (A1-DLT*A2) 

P2 = P1 

P3=SCALE*TAK (A1-ULT*A2) 

R Z = R 1 
p = Fl'**2 

IF(DLG,LT.O.) GO TO 1 1 

If(uLG.GT.18D. . ANO.DLG. LT. 180.) SO "O 41 
R1 = -P1*COS (OLG*A3) 

R2 = F2*C0S (DLG*A3) 

COMPUTE ACTIVE MAP MATRIX REGION 

V 1 = AN3-F 1- 1 . 

V2 = A.N3 + P2«-1 . 

IV1 = MAX1 (1. , VI) 

IV2 = HN1 (N1*1. ,V2) 

CALL SECOND (TX) 

WRTTF(6,500) 3,TX 

IF {TV1 . LF.IV2) GO TO 20 

IEXC=IV1 

I V1=MA.X0 (1 ,IV2-2) 

IV2 = MIN0 (N1,IEXC + 2) 

SO TO 20 

IF(ULG. GE. 270. .op. DLS. LE. DRi.) 3^ 21 

P2 = P3 

IF( (UI.S-180.) .LR". (180.-0LG) ) GO RO 30 
R2 = R2*COS (ULG*A.3) 

Gr TO 40 

IF ( (360 ,-ULG) . S’’’. DLG) GO TO 30 
P2=R2*CCS (ULGTA3) 

GO TO 40 

I F rULG. GE. 189 . . OP. DLG. I E. -1 80. ) SO TO ur 
IR (ULS . LT.90 . . AND. DLG. GT. -90. ) P’ =-”3 
IF (ADS (DLG) .GT. ULS) SO TO 13 
P1=P1*C0S(ULS*A3) 

SO TO 40 

P 1=P 1»ABS (CO S (DLG*A 3) ) 

GO TO HO 

COMFU'RE GRID BOX ARRAY ADDRESSES 

DC 100 I=IV1,IV2 
XH = 0. 

A P = R- ( AN3-I) »*2 
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fjofi OfJO orjfj ooo noo 


IF (A?) 100,1C3,1C4 

104 X!r=S0R?{AB) 

103 rH1=N3-XH-1 
I!I2=N3 + X?I*1 
CALI. FFADilS (4,P,S1 ,1) 
nr 200 .1=IH1,IH2 
FX = SC)PT( (AM3-.11 **2 + JAN3-I) *♦2) 
IF(BX.GT.BZ) GO TD 200 
11= (RLAT+2.*A4*ATAN (BX/SCALE) ) /XLAT+1 
IF(I1) 200,200,50 
59 AI=I 

IF(AI-AN3) 32,33,34 

32 AL0NG=ATA5i ( (AS3-J) /(AB3-I) ) ♦A4t130. 

GC TO 55 

33 ALONG=90. 


34 

55 


IF(J.LH. S3) AL0MG=270. 

GO TC 55 

AIOSG=ATAS ( (AH3-J) /{A S3-I) ) *A4 
TF (ALCSG.LT.O.) ALCSG= ALOSG + 3 60 , 
12= ( ALCMG-DLG)/XL0SG + 1 
TF(I2.GT.S) GO TO 200 
IF(I2) 200,200,42 


** PACK GRID BOX ARBAY ADDRESSES ISTO EACH HAP SATPIX SLEHENT 


p (.i)_=r2*i .E3+I1 
ON'^TSOE 


“2 , , 

200 COS'! 

CALL SPITHS (4,P,N1,I, 1,0) 

io(' crsTinnE 

CALL SECOND (TX) 
HRITE(6,500) 4,TX 
tt:=i 

IH1=N4 

15C CONTINUE 
IX=TZ+N4-1 

IF(IX.GT.H) IH1 = N4 + H-^IX 
♦♦ PFAD GRID FILS TAPE2 


1 5 


24 


43 


DO 15 J=1,IH1 

'JEAD(2,400) (DS (K, J) ,K=1 ,N) 
no 220 J=IV1,IV2 
CALL READHS (4,P,N1, J) 

DP 230 K=1,N1 
IF(P(K)) 230,230,24 
T2 = P fK) *1 .F-3 
I1=PfK -I2E1.E3 
IFII2) 230,230,43 

IF Ul. GT.IX.0R.I1 .lt. IZ) SO TO 230 
n = I1+H4-IX 


*• PROJECT PROPER GRID VALUES TO HAP HAT’RIX 


P(K)=DS (12,11) 

IF (DR (12,11) .EQ. 1. E5) P (K) =FLAG 
230 CONTINUE 

CALL WRITHS (4,P,N1,J,1,0) 

220 CONTINUE 
IZ=IZ+N4 

IF(IX.GE.H) GO TO 151 
GO TO 150 
151 CONTINUE 


CALL SECOND (TX) 

«PITS(6,500) 5,TX 

DO 160 1=1, N1 

CALL READHS (4,P,N1 , 1) 

DO 170 J=1,N1 

IF(P(J) .GT.100.) GO TO 145 
IFb(J)-FLAG) 170,145,170 


** FLAG NULL HAP FLEHFNTR 


145 P (J)=1 .E5 
170 CONTINUE 

CALL WRITHS (4 ,P,N1 ,1, 1 ,0) 

160 CONTINUE 

CALL SECOND (TX) 

HRITE(6,50C) 6,”^X 

** INVOKE COLORS ROUTINE 

CALL COLORS (N1,N1,P,IFILL, INDEX, 4, I”) 

CALL CLCSHS(4) 

CALL RECCND(TX) 
tfciTiJfe, 500) 7,TX 
400 'TRH AT (1CF1 1 . 4) 

5U0 FOPN AT (9H0STEP NC.,I3,27H IN POU'^INE P''LAP COHP LET E D , FR . 2 , 
*36H R'^CONDS AFTER START OF HA.IN PRi^oRAH) 

601 FCPHAT(2I4,4F1C.4) 

RFT'UPN 

END 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c +* 
c 


c 

c ** 
c 

33 


c 

c ♦♦ 
c 

1 

2 

3 

1 0 


c 

c »* 
c 


SUBRaUTI NE CCLORS (N 1 . N2 . V . IF, INDEX. lUN IT, IR ) 


THIS ROUTINE PSRFCPMS A USER CONTROLLED CONVERSION OF MAP PROJECTED SAR<? 
data into integer values RANGING FRCM 0 TO 21 IN UNIT STEPS. THESE INTEGER 

'"''■-SE COLOR SCALE AND hIvE A SELECTABLE LINEAR. 
exponential RELATIONSHIP TO THE SASS DATA. USER SUPPLIES 
Tn^A MA V LOVtER DATA VALUES AND NUMBER OF COLORS NEEDED UP 

CARDS^*^*^^^ COLORS. REFER TO MAIN PROGRAM FOR DESCRIPTION OF INPUT 


INTEGER UPOGUND.CNROUNO 
DIMENSION V (N2 >, IF( N2 ), INDEX { IR ) 

PRINT 39SS , 8 

COLOR SCALING AND DATA 8CUNCS INPUTTED 

REACIS.*) ISC ALE, ICOLOR.IBOUND.SIGHI.SI GLO 

ASSIGN 23 TC UP8GUND 

ASSIGN 23 TO DNDCUND 

I F ( I BCUND .NE ,2 I GO TO 33 

ASSIGN 21 TC UPDOUNO 

ASSIGN 22 TC DNBQUNO 

SET MAX NO, OF COLORS 

I F ( ICCLCP ,GT. 1 8 > IC0L0R=18 

IF( ICOLCR.GT.20) ICOLOP=20 
XCLR=I COLOP 
?CLP=XCLP+l .0 
I8NDMl=iaCUND-l 
I TOR-= ( I COLOF+I 8NDMI *2 ) 

COMPUTE DATA INCREMENTING VALUE 

IF( ISCALE-2 I 1,2,3 
DSTG=( SIGHI-SIGLOI/XCLR 
ASSIGN A TO JAIL 
GO TC 10 

DSIG=XCLR/AL0G(SIGHI-SIGL0+1 .) 

ASSIGN 5 TO JAIL 
GO TC 10 

OSIG=I SIGHI-SIGLO) /ALOGI XCLP+1. ) 

ASSIGN 6 TO JAIL 

DO 200 I =I , M 

CALL PEAOMS ( lUNIT , Y ,N2 . I > 

DO 250 J=1 , N2 

I F( Y( J ) . GE. S IGHI ) GO TO UPB CUND , ( 2 I . 2 3 ) 

IF! Y ( J ) .LE.S I GLO ) GO TO ONB CUN D , ( 22 , 23 ) 

COMPUTE INTEGER COLOR VALUE 


C 

C 

C 

A 

C 

C * A 
C 

5 

C 

c ** 
c 

6 

C 

C *♦ 

C 

21 


C 

C 

c 

22 

C 

C ** 

c 

23 

250 

C 

c ** 
c 


GO TO JAIL. (A, 5. 6) 

LINEAR CONVERSION 

I X-= < INT ( ( Y( J l-S IGLC )/OS IG) + I30UND J 
GO TO 250 

LOGARITHMIC CONVERSION 

I X= INTI DSIG* ALOGI Y( J )-S IGLOTl . ) )I- 18 CUND 
GC TC 250 

FXPON'^NTIAL conversion 

IX=(tNT(EXP((Y(J)-SI3Ln)/0SIG))+I8NDMI) 
GO TO 250 

FLAG HI DATA 

I X= ITOP 

I F ( Y ( J ) , EO. I . E 5 ) GO TO 23 
GC TC 250 

FLAG LO DATA 

I X= 1 

GO TO 250 
flag NULL DATA 
I X-0 

IF I J ) = I X 

INVOKE PACK ROUTINE 
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uuu 


CALL PACKEM< IF.N2.N3) 

♦* WRITE FALSE COLOR VALUES TO FILE TAPES 

buffer out <a.U ( IFI l » ,IF( N3) ) 
IF(UNIT(S>) 200. 301. 301 
200 CONTINUE 
RETURN 
301 WRITEI6.400) 

RETURN 

400 FORM AT ( I SHOBUFFER OUT ERROR. //IX) 

3999 FORMAT! 16H0ENTER CARO NO. .11) 

END 



CKJO r30«a OWfJ(3<Jf3(J 


SIinRf'IJTINE PACKBr" {IF,S'2,?:?) 




THIS RCUTINF PACKS TH? IJITBGBP AP^A^ I'='('!2) SnCM THA” BACH fiO 3IT WORD 
IS '’BDHCFD TO 8 PITS AND THESE 8 BITS APB dacKOD BACK INTO ARRAY IF. 
THE EFFECTIVE SIZE CF APPAY IF IS FFDUCFD TO M?=N2/7.5 


INTEGEF IF(N2),IB(B) 

TX=0 

** SET NO. OF 60 BIT WCR DS M EF DFD 


X = fI2/7. 5 + .9B9 
H3 = X 

«iu=»n/ 2 +i 

DO in 1 = 1, N4 
DO .-iO J=1,8 
30 IB(J)=G 
iy=IXAl 
K=(I-1)*15 

** PACK FIRST 7 1/2 (8 BIT) WORDS T"'"s 60 BI'^ WC’D 

DC 68 r,= 1 ,7 
KL=KH 

IF (KI.ST.N2) GO TO 22 
N= (P-I.) *8-4 
IP a) =IF (KL) .AND. 3773 
6? IB (t) =SHTFT (TB (I) ,N) 

IF ( (KL+1) .GT.H2) GC TO 22 
IB >8) =IF (K+8) . AMD. 360B 
IP(P) =SHIFT (TP (R) ,-4) 

IF (IX) =IB (8) .OP. IB (7) . OR. IB (6) . 0". TP (51 .DR . ID (4) .O'- . 13 (3) .CP. 
*IR (2) .OE.fB (1) 

DC 40 J=1,8 
40 IP(J)=C 
IX=TX+1 

^ + R 

iB(8) =IF(K) .AND. 17B 
IB (8) =SHIFT (13 (8) ,56) 

** PACK NEXT 7 1/2 (8 BIT) Wnpof; ^'iTr so BIT WnRR 

DC 75 L=1,7 
KI=K+L 

IF (Ki.GT.N2) GO TO 22 
N= (7-i) *8 

IP (I) =IF (KL) .AND. 3773 
75 IB (L) =SHIFT (IB (LI ,H)_ 

22 IF(TX)=IB(8).0P.IB(7).0P.IP(6) .OP. TP (5) .DP.IR(4).CP.IB(3).OR. 

*IB (2) .0?. TB (1) 

10 CONTIAMIE 
RF"’BRM 
END 
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PROGRAM ANNCT( INPUT, OUTPUT .T APE5= INPUT, TAPE6=OU TPUT » 


♦♦♦ THIS SOFTWARE WAS CREATED ON MAY 16, IS79 BY CUAYTON B. JACKSON WHILE 
WW* UNDER THE EMPLOYMENT OF THE RESEARCH TRIANGLE INSTITUTE CF N.C,. IT 
**♦ WAS DESIGNED TO SATISFY IN PART THE CONTRACT REQUIREMENTS UNDER NASA 
CONTRACT NASI-1E338. 


THIS PROGRAM CREATES AN ANNOTATION PLOT FILE CONTAINING DISPLAY 
LABELING AND A MAP LEGION, THE PLOT FILE IS A SPECIAL RASTER FILE 
AND IS COMPATIBLE WITH THE DISPLAY FILES GENERATED EY FALSCLR AND 
MAPGEN. 

THE INPUT CAROS ARE AS FOLLOWS: 

CARO I, ENTER A 10 CHARACTER MAP CODE I FOBMAT(AIO) I 

CARD 2, ENTER ICOLOR , I BOUN D , I TYPE . S I GH I , SI GL0 , AND XSCALE WHERE: 

I COLOR = NUMBER OF DISPLAY COLORS 

leOUNC = 0 - UPPER AND LOWER BOUND DATA NOT DISPLAYED 

1 - UPPER AND LOWER BOUND CATA DISPLAYED 

ITYPE = 1 - LINEAR COLOR SCALE 

2 - LOGARITHMIC COLOR SCALE 

3 - EXPONENTIAL COLOR SCALE 

SIGHI = UPPER BOUND CN DATA VALUES 

SIGLO = LOWER BOUND CN DATA VALUES 

XSCALE = PLOT DIMENSIONS IN INCHES 
*** FORMAT IS LIST DIRECTED ♦♦♦ 


DIMENSION XNUM(21 ),MC{ 1 ) 

CALL PSEUDO 

♦♦ ENTER MAP CCDE 

READ(5, 100) MC 
100 FCRMAT(AIO) 

** ENTER CONTROL VALUES 

READ ( 5,* ) ICOLOR , I BOON D , I TYPE , S I GH I , S I GLO , XSC ALE 
INX= ICCUQRHl 

I*^! IBOUNC,EQ.l ) INX = INX-2 

** COMPUTE COLOR BCUNCARY VALUES 

JNX=INX-1 
X NUM ( I NX )=S IGHI 
IF(ITYPE-2) 5,6,7 

5 OSIG=(SIGHI-SIGLO)/JNX 
GO TC a 

6 0S1G=JNX/AL0G( SIGHI-SIGLOHl.) 

. GO TO S 

7 OSI G=< SIGHI-SIGLO)/ALOG(FLOAT( INX ) ) 

3 CONTINUE 

DO 30 J=l.JNX 
GO TO ( 1 ,2, 3 ) . I TYPE 

1 '-INX-J 

XNUM( I ) =SIGHI-J*OSI G 
GO TO 30 

2 XNUM( J)i=SIGLO+EXP( ( J-1 )/OSIG)-l. 

GO TO 30 

3 XNUMI J)=SIGLO+OSIG*ALOG(FLOAT( J) ) 

30 CONTINUE 

SCALE'sXSCALE/S 12. 

CALL CALPLTISCALE.O ,5) 

XRECTp ( 256.-IC0LQRT1 2. )*SCALE 

YPECT= 1 5.TSCALE 

RW=24.*SCALE 

RH=1P.TSCALE 

** PLOT COLOR rectangles 

DO 10 I=l,ICOLCP 

CALL CALPLTl I , AHMASK,4 ) 

X=XKECT + (I-n*2A. ‘SCALE 

CALL RECTIX, YRECT.RH, RW ,0 .0 ,3 ) 

1 0 CONTINUE 

CSIZE-=7. ‘SCALE 
CALL CALPLT(63.0.4) 
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“4 


c 

c ♦* 
c 

15 


30 


40 


55 


C 

C ♦* 

c 


65 

75 

35 

95 


XNCT=XPECTMO. +SCALE 
VN0T=VRECT + 4.*SCALE ••• RH 
X0FSET=3.*SCALE 
IUP= I 

I S7=(ICOLQR+l)/2 

JF( laOUND.EQ.O) GO TO 15 

IUP=0 

I ST=< I ST/2) *2 


PLOT COLOR BOUNCAPY VALUES 


DC 30 I=l , 1ST 
IDX=( I-l (42+ 14- 1 UP 
X=XNOT+( I -1 ) 448 ♦♦SCALE 
X 1 = X 

I F( ABS( XNUM( IDX)).LT,10#) Xl=XlfXOFSET 

IF(XNUM< lOX ) .GE#0 . ) X1=X1+XCFSET 

CALL NUMBER (XI . YNOT , CS T 7E. X NUM < lOX) ,0 ♦ ♦ I ) 

CONTINUE 

ISP=INX-IST 

YNOT 1= YR ECT- ♦♦ 4SCALE-CSI 3E 
XN0T1=XRECT-14,*SCAL5 

IF( I BOUND. EG. 1 ) XNOT I =XNOT I +48 . *SC AL E 

DC 40 1=1, ISP 

IDX=142-IUP 

X = XNOT l+(I-l)*48.»SCALE 
X1=X 

I F ( A8S( XNUM< IDX)),LT,10.) Xl=X 1 +XCFSET 
I F ( XNUM ( I DX ) .GE.O. ) X 1=X1 +XCFSET 
CALL NUMBER (XI , YNOTl . CS IZE. XNUM< IDX ) , 0 . . 1 ) 
CONTINUE 

IF( I8OUND.EQ.0 ) GO TO 55 
XD I S=9 . *SC ALE 


Y=YN0T1 +3.*SCALE 

CALL PAPPO* (XRECT+XCIS.Y.XRECT-XDIS.Y, 1,0.0) 
IF( ISP. GE. 1ST) Y=vnOT+ 3.*SCALE 
X=XRECT+ICOLOR*24.*SCALE 
CALL PARPOW ( X-XDl S . Y , X+XDIS ,Y, 1 ,0 .0 ) 

CCNTINUE 

X=40.4SCALE 

IF( ITYPE.EQ . 1 ) X=56.4SCALE 
Y=55.*SCALE 


PLOT DISPLAY LABELING 


RADAR CROSS SECT I ON , 0 . , 30 ) 

o5«'5*85 


CALL NOTATE (999 . ,999. ,CSI ZE.21H 
GO TO 95 

CALL NCTATE(999.,999.,CSIZE,26H 
GO TO 95 

CALL NOT ATE ( 999 . , 999 . ,C S IZE ,26H 
CCNTINUE 

CALL N0TATE(999. ,999. .CSIZE ,5H 
CALL NCTATE(999.. '^99.. CSIZE, MC. 
CALL CALPLT(0. ,0. ,999) 

STOP 

END 


(LINEAR SCALE <DB)),0,,21> 

(LOGARITHMIC SCALE (DO)), 0., 261 

(EXPONENTIAL SCALE (CB)),0.,26) 

MC=,0.,5) 

., 10 ) 
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PROGRAM MAPGEN( INPUT. OUTPUT .T APE5=I NPUT t TAPE6=0UTPUT. TAPR3 . TAPES ) 


*** THIS SOFTWARE WAS CREATED ON JUNE 5. 1979 3V CLATTON 8. JACKSON WHILE ♦** 

UNDER THE EMPLOYMENT OF THE RESEARCH TRIANGLE INSTITUTE OF N.C.. IT *** 
«** WAS DESIGNED TO SATISFY IN PART THE CONTRACT REQUIREMENTS UNDER NASA 

CONTRACT NASI-1533S. AW* 


THIS PROGRAM GENERATES EITHER A MERCATOR OR POLAR STEREOGRAPH IC MAP PLOT 
WHICH CONTAINS LATITUDE AND LONGITUDE LINES. KEY MERIDEANS AND POLES. AND 
SUPERIMPOSED LAND/WATER BOUNDARIES. THE PLOT OFFSETS AND SCALE FACTORS 
ARE CALCULATED SUCH THAT THE PLOT FILE IS DIRECTLY COMPATIBLE WITH THE 
DISPLAY file created BY PROGRAM FALSCLR. 


THE CNE INPUT CARO IS DESCRIBED AS FOLLOWS: 

ENTER ZSCALE. MAPTYPE, ULT , OLT , ULG. AND OLG WHERE? 

ZSCALE = UNSCALED PLOT SIZE (IMAGE SIZE IN INCHES) 

MAPTYPE = I - MERCATOR MAP 

2 - POLAR STEREOGRAPHIC 

ULT = NORTHERN MOST LATITUDE (MAPTYPE-l) 

= OUTERMOST LATITUDE <MApTYPE=2) 

DLT = SOUTHERN MOST LATITUDE (=0 IF MAPTYPE*2) 

ULG = EASTERN MOST LONGITUDE (=0 IF MAPTYPE=2) 

DLG = EASTERN LEAST LONGITUDE (=0 IF MAPTYPE=2) 


♦♦ INPUT CARO USES LIST DIRECTED FORMAT ** 


DATA A1 . A 2. A3/. 785398 I 6 35. 8. 7266 A62609E-3 ,57.29 5779 31 / 

*♦ ENTER CCNTRCL VALUES 

READ(5, ♦) ZSCALE.MApTYP.ULT.DLT.ULG.OLG 
IF( MAPTYP.GT. I ) GO TO AO 

♦* FIND CENTER LATITUDE 

CL AT=2.»A3WATAN( SOFT ( TAN( A 1 HULTWA 2 ) *TAN ( A I HDLTW A2) 1 ) -90. 

** SET ANGULAR LIMITS 

P1=ABS(ULG-DLGM'0.5 

P2-PI 

P3 = 4BS ( CL AT-DLT ) 

P4=ABS(ULT-CLAT) 

PG=ULG-P1 

♦* PLOT MERCATOR 

CALL SUPMAP(9,CLAT,PG,0 ..PI .P2.P3.PA.-4. lO.-O.O. lEP) 

IF( lER.NE.O) GO TO 50 

♦* COMPUTE MAP IMAGE PIXEL SIZE 

0IJ=4.*A2*P2 

OV=ALaG(TAN{AlHP4*A2))+ALaG(TAN(Al+P3*A2)) 

IF(OU.GT.OV) GO TO 41 
M = 400 

N=400.*DL/OV + 0.5 
GO TO 61 
41 N=400 

M=400.*OV/DU + 0.5 

61 WRITE(6,5000) ULT . OLT , ULG , DL G. N . M 

GO TO 60 

** SET CENTER LATITUDE TC NORTH OR SOUTH POLE 
40 CLAT=89 .99999 

I F (OLT. LT. 0 . . AND.ULT.LE. 0 . ) CLAT=-89 .99999 
P t=ABS(CLAT-UL T ) 

♦* PLOT POLAR STEREOGR APH I C 

CALL 3URMAR<l,CLAT,0.,0.,Pl.Pl.Pl,Pl.-4.19.-0,0.!EP) 
IF(IER.NE.O) GO TO 50 
N = 400 
M = 400 
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oon nno 


6C 

** 

♦ * 

50 

1 00 

5000 

5001 


WPIT5( 6. 500 n N.M 
CALL FRAME 


COMPUTE OFFSETS ANO SCALE FACTORS 
SCALE=ZSCALE/51E. 

X SC ALE = A MAXI <10. *N/A608 .. 1 0.^M/A608.1 
YSC ALE-XSCALE 

XOFF = 0 .5*(5ia-N )*SCALE-0, 5*XSCALE 
YOFF=l 00.*SCALE-0.5*YSCALE 


CALL PLOT TRANSLATbR 

CALL METR3T < XOFF , YOFF , X SC AL E , YSCA LE .SCALE ) 

STOP 

WRITECe.lOO) lER 

CALL SYSTEM(52. 12HSUPMAP ERRORJ 
STC? 

FnPMAT(5H0IER=,I5) 

FuFMATCSAHOA MERCATOR MA° PLOT WAS GENERATED IN THE 
♦//6X.5HNLAT=.F6.2.5X.5HSLAT=,F6.2.5X,7FEL0NG1=,F7.2 
*F7,2 , //6X, 1 OH^hF map IS.I4.15H PIXELS WIDE BY.I4.7H 

FORMAT (48H0A POLAR STEREOGP APHI C MAP WAS GENERATED 
♦15H PIXELS WIDE BY.I4.12H PIXELS LONG) 

END 
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.5X , 7HEL0NG2 
PIXELS. 

THAT IS. 14. 
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PRCGRAM MERGE( INPUT .OUTPUT ,T APE5= INPUT . TAPE6=aUTPUT • TAPES . 
♦ TAPEIS.TAPE44) 


*** THIS SOFTWARE WAS CREATED ON MAY 14. 1979 BY CLAYTON B. JACKSON WHILE ♦♦♦ 
♦♦♦ UNDER THE EMPLOYMENT OF THE RESEARCH TRIANGLE INSTITUTE OF N.C.. IT ♦♦♦ 
♦♦♦ WAS DESIGNED TO SATISFY IN PART THE CONTRACT REQUIREMENTS UNDER NASA 
*** CONTRACT NAS1-1S23S. 


THIS PROGRAM IS DESIGNED TO MERGE TWO RASTER DISPLAY FILES BY PERFORMING 
A BIT BY BIT OR ON EVERY COMMON BIT. ONE FILE MUST HAVE 512 RECORDS. WHILE 
THE OTHER FILE SIZE IS USER INPUTTED. THE RESULTING FILE IS A RASTER 
DISPLAY FILE WHICH IS THE SUPER IMPC S I T ION OF THE TWO INPUT FILES. 


THE ONE INPUT CARO IS AS FOLLOWS: 

ENTER N AND M WHERE I 

N = HORIZONTAL PIXEL SIZE OF FILE 2 INC. OF 8 BIT WORDSI 
M = VERTICAL PIXEL SIZE OF FILE 2 (NO. OF RECORDS) 

** LIST DIRECTED FORMAT ♦♦ 


INTEGER BUFFI < 100) ,BUFF2( 100) 

NWOS=100 
PRINT 100 

♦4 ENTER 2ND FILE SIZE 
REACIS.*) N.M 

IF ( N.EQ.512 .ANO.M .EO.S 1 2) GO TO 79 

♦* COMPUTE BIT OFFSETS 

IXOFST=(5l2-N)/2 - 4 
I Y0FST=412-M 
I XD I TS= IX0FST48 
I ZB ITS=IX8ITS-( I XB ITS/60) 460 
I YBITS=60-IZBITS 
I XWRO=I XBITS/6041 

44 SET MASKING VARIABLES 

XMASK=.NOT.MASK(IYBITS) 
ymask=.nct. mask ( I zaiTs) 

ZMASK=MASK( IZBITS) 

IFIS12-N-IXCFST) 99.50.50 
50 IF( I YOFST.LT. 0) GO TO 199 
1F( lYCFST.EQ.O) GO TO 15 

44 PERFORM Y OFFSET 

DO 10 I-=l.IYOFST 

BUFFER IN (15.1) ( RUFF 2( 1 ). BUFF 2( NWOS ) ) 

IFIUNITIIS)) 20.299.299 
20 NWDS = LEr4GTH< 15) 

SUFFER OUT (44.1) ( SUFF2 ( 1 ) . BUFF2 ( N WD S ) ) 

IFIUNIT<44)) 10.399.399 
10 CONTINUE 
15 IDUMMY=0 

44 PERFORM MERGING - FILE SIZES UNEQUAL 
DO 150 1=1. M 

BUFFER IN (15.1) (BUFF2( I ) . BUFF2( NWDS ) ) 
IF(UNIT(15)) 60.299.299 
60 NWDS=LENGTH< 15 ) 

BUFFER IN (3.1) ( BUFF 1 ( 1 ). BUFFI ( NW OS ) ) 

IF(UNIT(8)) TO. 299. 299 
70 MWDS=LENGTH(8) -1 

I DUMMY= SHIFT (BUFF 1(1) , -IZBITS). A. YM ASK 
BUFF2< IXWRD )= IDUMMY.OR.aUFF2< IXWRD) 

DO 140 J=l.MWOS 

44 PERFORM X OFFSET 

ITEMP l = eUFF I ( J ) . A.XMASK 
ITRMPl=SHIFT(ITEMPl.IYnlTS).4.ZMASK 
I TEMP2=SHIFT( BUFF 1 ( J4 1 ) I Z B I T S ) . A . YM AS K 
IDUM:MY=ITEMP1 .OR.ITEMP2 

44 PERFORM THE BIT BY BIT OR 
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BUFF 2 ( IXWRO-*- J) =eUFF2( IX \«R0 + J ) . OK. I DUMMY 
140 CONTINUE 
J ■— C S I 

IF{ ITBITS.EQ.O) go to 144 
ITEMPlaBUFFl ( J) .A.XM4SK 
I DU MMY= SHIFT! I TEMP I ,I V01TS) .A.ZMASK 
BUFF2( tXV»RD-l-J ) = BUFF2( I X V>PD+ J ) . OR ♦ I D UM M V 
144 BUFFER OUT (44.1) ( BUFF2 ( 1 ) , 3UFF2 < NWDS ) ) 

IF(UNIT<44)) 150,359,359 
150 CONTINUE 

*♦ COMPLETE OUTPUT FILE 

34 BUFFER IN (15,1) ( BUFF2( 1 ) ,8UFF2( NWDS ) ) 

IF(UNIT(15)) 130,230,299 
130 BUFFER OUT (44.1) ( EUFF2( I ) , BUFF2 ( NWDS ) ) 

IF(UNIT(44)) 34,399,399 
230 ENOFILE 44 
ENDFILE 44 
REWIND 44 
REWIND 15 
REWIND 9 
STOP 

♦♦ PERFORM MERGING - FILE SIZES EQUAL 
79 DC 250 1=1,512 

BUFFER IN (15,1) ( BUFF2( 1 ) , BUFF2( NWDS ) ) 

IF(UNIT(IS)) 330,299,299 
330 NWOS=LENGTH( 15) 

BUFFER IN (9.1) (BUFFI ( 1 ) .BUFFI (NWDS) ) 

IF(UMT(8)) 340.299,299 
340 DO 260 J=l,NWOS 

260 BUFF2(J)=BUFP2( J ) .OR .BUFFI ( J ) 

BUFFER OUT (44,1) ( BUFF 2( 1 ) . BUFF2 ( N WO S ) ) 

IF(UNIT(44)) 250,399,399 
250 CCNTINUE 
GO TO 230 

99 WRITE(6,200> IXCFST 
STOP 

199 WRITF(6,300) M 
STOP 

299 WPITE(6.400) 

STCR 

399 WRITF(6,500) 

STOP 

100 FORMAT (29H0ENTER N , M TFE MAP PIXEL SIZE) 

200 FORMAT! IHO, 35HX-OFFSET IS SUCH THAT DATA WILL NCT, 

♦ 18H FIT IN PLOT FRAME . 5X, 7H I XOFST= . I 4 , 7H PIXELS) 

300 FORMAT(31HOY-OFFSET IS LESS THAN ZERO — M=,I4.12H IS TO LARGE) 

400 FORM AT ( 25H04* ERROR ON BUFFER IN ** ) 

500 FCRMAT(26H0** ERROR ON BUFFER OUT **) 

END 
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